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ABSTRACT 

This work is an attempt at determining the optimum 
neutral axis location in the midship section of a vessel, 
using as a model for the calculations a box girder, subject 
to several different conditions of sagging and hogging 
bending moments^ several values of design stress are also 
considered and in the end what happens when lateral hydrosta- 
tic load is considered along with different values assigned 
to the design stress in tension and in compression is also 
analysed. 
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a 


- distance between transverse stiffening members 


AL 


- length of box girder or vessel 


B 


- beam 


BOT 


- thickness of the plating on the bottom 


BM 


- bending moment 


D 


- depth 


DIBOT 


- distance of the bottom from the neutral axis 


DITOP 


- distance of the top from the neutral axis 


DSC 


- design stress in compression 


DST 


- design stress in tension; also design stress without 
discriminat ion when DSC=DST 


H 


- head of sea water ( lateral load ) 


HB 


- head of sea water on the bottom 


HM 


- hogging bending moment 


HT 


- head of sea water on the top 


I 


- moment of inertia of the cross section with respect 
to the neutral axis ( GMI in the program ) 


LBS 


- longitudinal bending stress 


S 


- spacing of longitudinal stiffening members 


SBCRIT 


- critical stress on bottom 


SBH 


- stress on bottom due to hogging 


SBS 


- stress on bottom due to sagging 


SM 


- sagging bending moment 


STCRIT 


- critical stress on top 


STH 


- stress on top due to hogging 


STR 


- stress 



STR1 


- 


primary stress = LBS 


STR2 


- 


secondary stress 


STR3 


- 


tertiary stress 


STR3B 


- 


tertiary stress on bottom 


STR3T 


- 


tertiary stress on top 


STRCR 


- 


critical stress 


STS 


- 


stress on top due to sagging 


t 


- 


thickness of theplate 


T 


- 


draft 


TOP 


- 


thickness of the plating on the top 


TS 


- 


thickness of the side plating 


TSTH 


- 


total stress on top in hogging 


TSTS 


- 


total stress on top in sagging 


TSBH 


- 


total stress on bottom in hogging 


TSBS 


— 


total stress on bottom in sagging 



INTRODUCTION 



7 



It is frequently found in the literature related with 
the strenght of ships such statements as: 

"In general the neutral axis is expected to be at 40$ of 
the depth (D) above the base line ... a preliminary estimate of 
the thickness required for the bottom plating (BOT) can be 
obtained from 



0.4D 3/ 40300 

1.875( ) 2240 \ AL = 

0.6D V 1.75 

1 S 

1 + ( ) 

950 BOT 



where the Montgomerie's expression is being considered to 
determine the critical stress (SBCRIT psi) and the short edges 




a- distance between transverse stiffening members (inches) 
1.875 accounts for a factor of safety 

Tobin's expression was considered in the above formula 
to give the maximum permissible bending stress (LBS, psi) 

J / 



LBS= 2240 



AL 
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We see thus that this procedure allows for an estimate of 
the bottom thickness to be obtained after knowing the length and 
assuming a longitudinal frame spacing. What may be questioned 
here is the reason why one assumes the neutral axis 0.4D above 
the base line and also the reasons behind the choice of a frame 
spacing. 

In the following a procedure to estimate the scantlings 
of longitudinal members is described; this procedure can also 
be found frequently in the literature dealing with the subject. 

. . . "For a particular set of principal dimensions , general 
arrangement and frame spacing, plate panel thicknesses for one 
after another of the amidships longitudinally continuous elements 
are selected for an assumed position of the neutral axis at 40^ 
of the depth above the keel. The composite midship section assem- 
blage is thentested for the adequacy of its section modulus when 
ranged alongside some standard value, such as given by the 
Load Line Regulations or some function of it. Thus the modulus 
is judged sufficient to limit ship bending stresses to the value 
previously set and incorporated into the design equations. 

For all but the deepest ships the modulus will undoubtedly 
need to be increased and a small arbitrary increment is added 
to the upper deck or bottom whichever happens still to be farther 
from the true neutral axis location found after the first design 
cycle. The shift of neutral axis resulting from the scantling 
increment requires a compensating increment in the opposite flan- 
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ge with the obvious outcome that the neutral axis is shifted in 
reverse with attendant effect on the scantling chosen first; a 
cyclic process is then begun. More trials will be necessary 
before the match of available to required section modulus is 
achieved. " 

Prom this transcript, we may again note the initial 
statement "...for a particular set of principal dimensions, 
general arrangement and frame spacing. . . " and question whether 
or not the frame spacing should be chosen as an input or rather 
be left as a variable and be obtained as a result of a structural 
design procedure. It is here aknowledged though that in these 
procedures an acceptable solution can be worked out and each 
new cycle is an improvement over the precedent ones. What may 
be questionable though is whether the first cycle, due to the 
assumptions made at its beginning, has or not a strong influence 
in the outcome leading to a solution among the various possible 
ones which in the end does not lead to"the optimum design" 
whatever its definition may be. 

It seems that this trial procedure can be looked at as a 
"weakness" of the structural design procedure for the determina- 
tion of scantlings and it is perhaps reasonable to make some effort 
in order to work out a process that could give the scantlings 
in a direct manner and be at the same time in agreement with 
what can be considered "the optimum design", (usually based on 
economic considerations ) . 
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Nowadays much concern is given to economics in ship design 
and naturally the strength calculations of a ship are also made 
with economical optimization in mind. With this respect, we can 
realize that in ships structures there is a desirability to mi- 
nimize structural weight and some define structural efficiency 
as being the least weight structure capable of resisting the 
loads for which they are designed. 

It should be clear that the cost of construction and of 
maintenance and repair are inter related factors together with 
least weight, to be considered when bringing economics into 
ship structural design. Also other matters such as accessibility, 
weldability and availability of material must be considered as 
well . 

In the following sections an attempt is made to obtain 
without any interactions the scantlings for the plating of a 
box girder. It is hoped that through this simple model some 
helpful conclusions can be made in view of their application 
to a vessel's structure. Each of the different situations ana- 
lized will be exhaustively explained in the following sections. 
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PRELIMINARY CALCULATIONS 

The model used in this section and also throughout this 
work with only very small additional considerations is a box 
girder which cross section is pictured in figure 1. 




As will be seen in the calculations of the following 
sections, the dimensional parameters used are of the order of 
magnitude of those considered in ship design. 

This model was chosen due to its simplicity for the 
performance of the following calculations as can be understood 
immediately from its geometrically symmetrical cross section. 

The objective of this section is to determine the neces- 
sary formulation to be used in future sections. 

The formulation dealt with in this section is related 
only to the particular model chosen and the geometry of its 



cross section. 
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£(A*d) 

DIBOT= *12 

Z A 

Z(A*d) 

DIBOT 1 

Z A 

Z A=B*T0P+2*TS*D+B*B0T=2*D*TS+B* ( TOP+BOT ) 

DIBOT =Neutral axis location above bottom (in) 

DIBOT' = " " " " " (ft) 

■^BOT =Moment of inertia with respect to the bottom 

I _ M n H n n n n neutral axis 

Z(A*d)= ( D/1 2 ) * ( B*TOP+D*TS ) 

I=Z(A*d*d)fZ I-( Z(A*d)) 2 /( Z A) 

DI BOT=D* ( B*TOP+D*TS ) / ( 2*D*TS+B* ( TOP+BOT ) ) 

This formula will give DIBOT in inches if D and B are 
expressed in inches; it will give DIBOT in feet in the event 
that D and B are expressed in feet. 

DITOP=D-DI BOT 

DITOP=D* ( B* BOT+D*TS ) /{2*D*TS+B ( TOP+BOT ) ) 

DITOP in units of inches or feet as explained for DIBOT 
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The following formula will give I in units of (in^ft)^ 
for B and D given in inches. 

( D/12 ) 2 * ( 3*B*B*T0P*B0T+( D*TS ) 2 +2*B*D*TS* (TOP+BOT ) ) 

1 = 

3* ( 2*D*TS+B* ( TOP+BOT ) ) 

Next formula gives I in units of (in*ft) for B and D 
given in feet. 

Z 2 

4*D*D* ( 3*B*T0P*B0T+ ( D*TS ) +2*B*D*TS* ( TOP+BOT ) ) 

1 = — — 

2*D*TS+B* ( TOP+BOT ) 

These expressions will be used in most of the following 
sections with B and D expressed in feet. 

NOTE: TS which is defined as the thickness of the 
side plating is considered throughout this work to be uniform 
and will always be equal to the larger of the thicknesses obtained 
for the top and for the bottom plating. 
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FORMULAS USED IN THE ANALYSIS 



BEAM EQUATION 



Assuming pure bending, the stresses in the cross section 
of the box girder are given by 



BM*y 

STR= 

I 



BM 




neutral 

surface 

x 



Fig. 2 

where y is the distance from the neutral axis to the point 
where the stress is being determined. 

The neutral surface indicated in fig. 2 is composed of all 
the points which are not subjected to any stress in pure bending. 



BRYAN EQUATION - Plate under inplane load 



(fl) 2 *E*(t) 2 

STRCR=K* j 

12*(1-(P0I) )*(S) 

V/hen a/S>2 then K=4 

STRCR given in psi for t andS expressed in inches and E in 

psi . 



a 



S 



=3. 1416 
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PLATE UNDER LATERAL LOAD 
S P 1 

STR3=0. 5*K*P*H*( ) * 

' t 144 

STR3 - tertiary stress 

- density of sea water (64 lb/ft-^) 

H - head of sea water (ft) 

a - long dimension of plate (in) 

S - short dimension of plate (in) = distance between 
longitudinal stiffeners 

t - plate thickness (in) 

The stresses are more severe at points A and B but we are 
particularly interested in the stresses at point B because there 
they will combine with the stresses resulting from bending also 
called primary stresses. 

Note that for a panel under lateral load, the direction of 
the stresses is, at the edges, perpendicular to them, so at B 
they are directed "longitudinally". 

For a/S>2 we have K=0.685 

Replacing all the numbers in the above formula, we obtain 

s 2 

STR3=0. 152222*H*( ) 

t 



1? 
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The following stress definitions and figures 3 and 4 will 
help understanding the reasons behind the calculations made in 
many of the following sections. 

STRESS DEFINITIONS 

STR1 - SHIP BENDING STRESS or PRIMARY STRESS 

Since this stress is directly proportional to distance from 
the ship's neutral axis, it is substantially constant at any point 
on the cross section of the bottom or top plating when the vessel 
is in an upright attitude. For distinction, this stress may be 
termed a primary stress or an area stress in this location and 
under such a condition. V/hile the longitudinal direction is imp- 
lied in the foregoing comments, they are equally applicable 
transversely. 

STR2 - GIRDER BENDING STRESS or SECONDARY STRESS 

Girder bending stress arises from the reaction of a plating 
stiffener combination to a loading normal to the plating such as 
may be due to v/ater pressure. Following the usual concept, the 
plating acts as one flange of the composite girder and so is 
stressed. Due to shear lag, this secondary stress in the plating 
is a maximum in way of the stiffener and diminishes with distance 
from it. 

For the previous case of bottom plating and vessel upright, 
the transverse variation of fore and aft stress will be sensibly 
constant for some distance from the stiffener but will vary 
significantly through the plating thickness as the plate thick- 
ness is a substantial part of the distance to the neutral axis 
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of the plate stiffener combination. For convenience such a stress 
may be termed a line stress although for a wide spacing of stiffe- 
ners with consequently greater deviation from uniform distribution 
some inaccuracy of concept is introduced by using this term. 

For girders continuous over several supports, such as tran- 
sverse bulkheads, maximum bending moments occur at the supports 
and will readily be in the order of twice the field moment. 

STR3 - PLATE BENDING STRESS - TERTIARY STRESS 

Proceeding in turn to ever more elemental portions of the 
ship's structure concludes with the simple plate panel supported 
on its four edges. In such a plate under lateral or normal loads, 
bending stresses are created which are variable from point to 
point on the panel and through its thickness with maxima occurring 
at one or the other surface. Since any such stress is effective 
at a point only, it may be called a point stress or tertiary 
stress . 

Each elemental panel of ship's plating is generally one of 
a repeating pattern. Quite appart from any rotational restraint 
provided by the panel supports consisting of floors, frames, 
longitudinals, etc., there is the restraint afforded by the 
continuity of the plating and its loading. For pratical purposes 
then, a condition of edge clamping exists with maximum stresses 
being found at the midlength of the longer panel edges and 
secondary maxima at the midlength of the shorter edges. In both 
cases these stresses are normal to their respective edges. 



BENDING STRESSES IN SHIPS 
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STRESS DISTRIBUTION IN 
BOTTOM SHELL AT 
TRANSVERSE WEB FRAME 
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HULL BENDING STRESS 



TYPICAL BOTTOM STRUCTURE 
OF A 

LONGITUDINALLY FRAMED SHIP 
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STIFFENER BENDING STRESS 
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PLATE BENOING STRESS 




FIGURE 3 
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From the previous definitions and figure 3 we may summarize 
in figure 4 the effect of STR1 , STR2 and STR3. 




Fig. 4 

At point A 

STR2 is larger than at point B 
STR2 is acting in the longitudinal direction 
STR3 is larger than at point B 
STR3 is acting in the transverse direction 
STR1 and STR2 are here combined acting both in the 
longitudinal direction 



At point B 

STR 2 may be acting in the tansverse direction 
STR3 is acting in the longitudinal direction 
STR1 and STR3 are here combined acting both in the 



longitudinal direction. 
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Usually STR2 is attributed a constant value and can be 
taken into account by automatically decreasing the design stress 
by the value assigned to STR2. 
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THE BENDING MOMENT EQUATION 



The bending moment equation used is derived from 



knowing that 



and that 



and using 



we obtain 



BM = A*AL/35 ton*ft 

A = V/35 ton 

V = B#T*AL*C b 

C b = 0.75 

0.75 

BM = *B*T*AL*AL ton*ft 

35^35 



where 



AL = is the length in ft 
B = is the beam in ft 
T = is the draft in ft 

= is the block coefficient 



or multiplying by 2240 we obtain 



BM = 1 . 37143*B*T*AL*AL lb*ft 



Note : 

V is the displacement volume and ^ is the corresponding 



displacement weight. 
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MOMENT ACTING IN ONE DIRECTION 

NO RESTRICTION IMPOSED ON THE COMPRESSION STRESS 

Considering the box girder subjected to a bending moment 
acting in hogging BM, we will observe as a consequence that the 
top will be subjected to tension and the bottom to compression. 
These stresses are given by the already introduced formula 

BM*Y 

STR= 

I 

or in this case, being interested on the top and bottom, 

BM*DIBTOP 

STH= 

I 

BM*DIBOT 

SBH= 

I 

In this introductory study it was decided not to limit the 
value of the stress in compression. Only the value of the stress 
in tension, here STH will be requested to be equal to the value 
chosen for the design stress in tension (or smaller). 

SBH will be a consequence of the value of DST and also 
of the S value chosen. 

So, we will basically have 3 equations 

DST = STH+*-TOP Bending equation 



S SBH BOT 



Bending equation + Bryan formula 
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BM 



STH 




COMPRESSION -* 1 ► TENSION 

Pig. 5 



The three equations that can be obtained will allow one to 
determine the three unknowns which are BOT,SBH,TOP. 

By replacing in the equations already mentioned the expre- 
ssions derived previously for DIBOT, DITOP and I and solving for 
BOT,TOP and SBH we obtain after a lengthy but simple process the 
equations in their final format as we will see. 

Note here that in order to solve the three equations, • 
the thickness of the side plates TS was made equal to the larger 
TOP or BOT; since TOP and BOT were not known and had to be deter- 
mined after deciding on TS, the equations had to be solved twice. 
First the equations are solved assuming that TOP> BOT this 
simply implying that the neutral axis is going to be chosen alrea- 
dy closer to the top and TS = TOP. Second, the equations are 
solved assuming BOT > TOP thus making TS = BOT. 

As will be seen, only one solution is possible, the other 
being rejected by giving a result that does not agreee with the 
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assumption just mentioned. 

TOP , BOT , S (inches) 

B , D , DITOP , DIBOT , AL (ft) 

DST , STH , SBH , E (psi) 

BM (ton#ft) 

I (in*ft ) 2 

560*BM* ( B*BOT+D*TS ) 

DST> STH = 5 - 

D* ( 3*B*B*TOP*BOT+2*B*D*TS* ( TOP+BOT ) = ( D*TS ) ) 

560*BM* ( B*TOP+D*TS ) 

SBH 

Same denominator 




(S) 2 *(1-(P0I) 2 ) 

G = 0.30396* 

E 

DST B*BOT+D*TS 
SBH B*TOP+D*TS 



Making 

9*B*B*D 

A =- 

l40*BM 

3*B*D*D 

C = 

?0*BM 

3 *d*d*d 

F 

l40*BM 



we will obtain, by assuming first TOP> BOT 
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D*(B+D)*TOP 

DIBOT = 

(2*D+B)*T0P+B*B0T 

D*(B*BOT+D*TOP) 

DITOP 

Same denominator 

( 2*B*D+D*D)*TOP*TOP+( 2*B*D+3*B*B)*TOP*BOT 

I = 4*D*D* 

Same denominator 

B* ( BOT ) 3 

TOP — 

( B+D ) *G*DST-D* ( BOT ) ^ 

TOP* ( 1 2*D- ( C+P ) *DST*TOP ) 

BOT 

DST* ( A+C ) *TOP-12*B 

Before going to the equations obtained when BOT > TOP 
is assumed, we call the attention of the reader to the fact that 
the values of BOT and TOP have to be obtained from the two above 
equations. Substitution of one into the other would result in 

4 - 

an equation of degree higher than at least 5 degree, unsuitable 
for direct solution. As will be later exemplified, the solution 
is here obtained by assigning values to either one of BOT and TOP 
and through the two equations, we obtain the respective values 
for TOP and BOT. The ploting of these values BOT vs TOP will 
give one line for each of the two equations. The intersection of 
the two lines thus obtained allows the solution to be read from 
the plot. 

This solution will be acceptable if the intersection lies 
in the sector TOP> BOT otherwise the intersection must be 
rejected. 

Assuming nov; BOT>TOP we make TS = BOT and obtain: 



2b 



d*(b#top+d*bot) 

DIBOT 

(2*D+B)#B0T+B*T0P 

D* ( B+D ) *BOT 

DITOP 

Same denominator 

( 2*B*D+D*D ) * ( BOT ) 2 + ( 2*B*D+3*B*B ) *TOP*BOT 

I = 4*D*D* — 

Same denominator 

( B+D ) * ( BOT ) 3_D*G*DST*BOT 

TOP — 

B#G*DST 

12*(B+D)-(A+C)*DST*T0P 

BOT 

( C+F ) *DST 

In this form the equations can be used in a very simple 
computer program that will give the set of values TOP vs BOT 
suitable for ploting. At this stage we are hoping that if the 
intersection obtained in the situation where we assumed TOP> BOT 
is acceptable , the two lines resulting from the case BOT> TOP will 
not intersect or if they do that the intersection must be rejected 
due to the fact that they will violate the assumption from which 
they were derived 

As can be observed from the several cases run in the compu- 
ter, the case where two acceptable solutions exist never happe- 
ned and it can only happen if TOP = BOT where then both solutions 
are acceptable and coincide as in the case of one of the follo- 
wing runs. 

In the following pages the computer printouts are included. 
The first run is explained in detail to make the interested rea- 
der familiar with the problems encountered in determining the 
intersections and the respective thicknesses. 



29 



The other runs are here presented as well to serve as 
proof that a simple solution can always be obtained and also to 
give an idea of the range of thicknesses expected in each case. 

In all runs only a bending moment acting in hogging was 
considered. 

The program used is a direct application of the equations 
here derived to a particular set of values chosen as the input. 
From the results obtained, some conclusions were drawn that a- 
llowed for the improvements made in following sections. 

No further explanation on the program shall here be given 
because it is self explanatory in the comment statements 
included. 
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From these lists of thickness values printed by the com- 
puter, we can see that the acceptable ones are those marked 



sigi 



with a II sign. Plotting these corresponding pairs of values 
properly, as is next exemplified, allows for one to read from 
the intersection of the two lines, which falls within the range 
set by the assumption that originated it, the values of thick- 
nesses at the top and at the bottom which satisfy the model 



we started with. 



35 




We can see that the thickness values are easily read from 
the intersection obtained and they are acceptable because they 
fall within the range assumed to derive the formulas that ori- 
ginated them namely TOP> EOT. 

Next will be plotted the case BOT>TOP and as can be seen, 
again an intersection is obtained which in this case is not 
acceptable because it occurs within the range TOP> BOT contrary 
to the assumed BOT > TOP that originated it. 

With a little more experience, this can be seen immediatly 
from the computer printouts, thus avoiding loss of time doing 
a plot that will lead to a nonacceptable result. 
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Now having obtained the values for the thicknesses from 
the plots, we verify them through the use of the next program 
based on the formulation previously derived. As can be noted, 
we obtain an acceptable value for the stress at the top (since 
we are considering a hogging bending moment) which we said 
before in the derivation of the equations should not exceed the 
design stress in tension. In this case the design stress was 
DST=20000 psi and we see that the stress at the top STH=20037.5 
the 37.5 coming possibly from inaccuracy when reading the thi- 
cknesses from the plot, and computer roundoffs. 

Note that the stress at the bottom, since there was no 
limit placed on it except the Bryan formula to resist buckling 
caused by compressive inplane loading, is attaining a too large 
value which is unacceptable and unrealistic. Here again there 
should be no difference between SBH and SBCRIT and the difference 
that can be noted is also due to inaccuracy in reading the thi- 
cknesses from the plot and computer round-offs. 
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SUMMARY FOR THE PRECEDING CASE 



Length 




500 ft 


Beam 
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Draf t 


S 


B /3 


Depth 


= 


AL /9 


DST 


= 


20000 psi 


S 


= 


10 inches 


Material 


“ 


STEEL 


TOP 


= 


0.468 inches 


BOT 


= 


0.166 inches 


STH 


= 


20037.45 


SBH 


= 


30434. 89 


SBCRIT 


= 


29886.80 


DITOP 


= 


22.06 ft 


DIBOT 


=S 


33.50 ft 



We can already notice here that the stress in the bottom 
is attaining a too large value due to the fact that, since no 
limit was imposed on it, it is dependent only on the chosen 
value for S and the respective capacity to resist buckling pro- 
vided by S and BOT while at the same time stressing the top plate 
to its maximum allowed DST. 

In the following two cases only the spacing S between the 
longitudinal stiffening members will be allowed to change and 
the results obtained are simply stated as previously exemplified 
and were determined in a similar way. 

From the results obtained we can note the effect caused 



by the change of S. 
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v r 



»' N S r F H ACCESS ?07O 
C J T ’.ON BEGINS: 



END 



STH= ? '037 .113 
S6H = 2T37.113 
SPCPIT= 2' 5/5 • 3a g 
0 I TO p = 27*7778 
OlBOT= 27*7778 



f v: 


C OPEPATIf.G SYSTEM 


1 vEpSION 


1 REVISION 


P12 03/0A/74 GENERATED 05/06/7A 


,00 


HANDLING CHARGE 


*• . 35 


/ 


JOB 


.35 


F5 


LINES PPI-iTEO PR2 


4 1.25 


/ 


X LN 


• 1 1 


5a 


CAROS read 


£ 1.50 


/ 


< CO 


.08 


?0 


RLOTTER VrCTOPS 


£ .25 


/ 


1/710 




13 


*• CDEL 70 seconds 


425.0'/ 


/ 


UplJR 


•09 


0? 


*■ ODEL EP SECONDS 


412.57 


/ 


HOUR 


• 0v3 






TOTAL 


CHARGE t 


•63 



. $ 



FElP a 9/ j a 7 3 1. LOGGED OUT /5/0B/7A 15*.a6 



11.1ft left after 



21 LOGINS 



SUMMARY 



Length 




500 ft 


Beam 


= 


AL/8 


Draft 




B/3 


Depth 


= 


AL/9 


DST 


- 


20000 psi 


S 


- 


30 " 


Material 


- 


STEEL 


TOP 


= 


0.4125" 


BOT 


= 


0.4125" 


STH 


= 


20037.11 


SBH 


- 


20037.11 


SBCRIT 


- 


20505.43 


DITOP 


= 


27.78' 


DIBOT 


= 


27.78' 



Note that this is the case where two acceptable solutions 
exist although they end up being the same because TOP=BOT. 

Note that the neutral axis is located half way between 
top and bottom and the stress is the same at the top and at the 
bottom, as it was expected beforehand. 



STRENGHT Of SHIPS-ON the neutral axis location 





c 




*A L * IS THE L ENGHT 


if 






A L * 5 0 0 • 


occ 






R = A l / 8 • 


}? 1 8 


c 




*DPA* IS ThE draft 


Cl? 






DRa=B/3 • 








D=AL/9 . 


i ' 3 U 


c 




*POl‘ IG ThE' POISSON'S RATIO OF THE MATERIAL 


\'?3V 


c 




for steel 


-3:i 






P 0 I «= 0 • 3 


7?H 


r 




*E« IS TmE YOUNG'S MODULUS OF THE MATERIAL 


“ 38 


c 




FOR STEEL 


~ 3 8 






E = 3 O • * 1 •*' • * * 6 


r^c 


/"* 

L. 




* D S T * IS The DESIGN STRESS 


sr 






D S T = 2 v’ 0 ■ ? • 


7 5 4 


r* 




* 3 •« IS THE SPACING OF LONGITUDINAL STIFPENING MEMBERS 


754 






S= 1 20 • 


•7cC 






IS the BENDING moment 


75^ 






' N '’=(B , -DRA x AL*AL*'/.7b)/(35.“35*) 








D = T • 3 0 3 9 6 * S * S * ( l.»c>cl*P3I ) / E 


r s4 






A — B**c’*B*D/( i 4 / • ** P ^ ) 


"DC 








1'^ 






f = 3 . - D - D * D / ( \ 4 m ) 


M?C 


r 

w 




THICKNESS OR TCP ASSUMED LARGER THAN BOTTOM THICKNESS 


:icC 






R I Tr ( p , 1 ? ) 


’ t <* fc 




10 


r3P v AT ( ' i ' , ?sX , ' tqp ASSUMED LARGER THAN BOTTOM') 


'•AC 


r 




.TCP* IS ThE T'.-ICKNESS of THf TOP PLaTinG 


;;fc 


r 




* P C T * IS T.jF THICKNESS OF THE BOTTOM PLATING 


r f: 






a'JT i =*.i . 


•7S 




i i 


-Co' =««Ro”l*f?OTi*3CTl/( (B + D) * G * DST - D* BoT 1 * BoT 1 ) 


* * o c 






W R I " E t 5 » ?0 ) T o D 1 > R n T 1 


• * r ^ 




20 


F^R^AT ( ' ' i 2 5 X / •Tr.p 1 -IF9. 4 ,5X, • »0T1 s ' Fg . 4 ) 


? m !> 






p OT = n OT i - . ?5 


?lf 






I r ( n QT l . qT .2. ) GO TO 30 








GO ”0 11 


•p?c 




30 


T^P ?*(’ • 0? 


p3* 




31 


’-JlT? = *O p P" ( 1 ? • * D ~ ( C*F ) *DST * T 0P2 ) / ( DST * ( A + C ) *T0P2**12 • *1 


I?*' 7 






VElTE<5>40iT0 p 2'BDT2 


‘? ? 4 




40 


FORMA’ » ' ' . ? r X / ' TpPP= ' F9 • 4 / 5X> 'B0 T 2='F9.4) 


I* cCr 






* c p ? = t o p ? ■* •: . ? 2 








l~ ( TOP? ■ GT . 1 . ) GO TO 50 


;f i 






GO ’0 31 


I'pPr 




50 


CON" IMJE 


? r r 


c 




T’ J ICINESS n F BC T TO^ ASSUMED LARGER THAN TOP 


• p r r 






WRI*p t 5, 60 ) 


31? 




60 


FORMAT {' l ?5X, ' OOTTOM ASSUMED LARGER T u A N TOP') 


■?*e 






pO t 3 = 1 • 


•np 




61 


T0p3= M P + D 1 *pOT3*ROT3*9OT3-D*G*DST*0oT3 ) / ( B*G*DST ) 


*334 






RR! T F ' 5/ 70 ) T oP 3 / B0 T 3 


3P p 




70 


rORMATl ' ' ,2hX, 'TCR3='F9*4/5X, 'B0T3='Fg,4) 


*?Ee 






B 0 T 3 = B 0 T 3 * C *75 


’3*C 






IF ( B 0 T 3 • G T .2 . ) GO TO 80 


p 3^e 






GO TO 61 
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STRENGHT Op Sl-'IPS-ON TwE NEUTRAL AXIS LOCATION page 2 

-a02 80 TOpA=e.e? 



'40 A 


81 


nOT<f=(12.*(B + D) 


-( A+r)*DST*T0P4)/( ( C + F > * 


DST ) 














WRITE ( 5> 9C ) T()P4 


/ POT A 
















)t7? 


90 


FORMAT < ' ' , ?9X, 


' TCPA 


c*F9*4»5X 


, ' BOT 4 = ' F 9 


. 4 ) 










H9 A 




TOPAeTOPA+0. 02 


















:aA£, 




IR ( TOP A • GT . 1 . ) GO TO 


100 














jab* 




GO TO 81 


















jaPc 


100 


CONT INUE 


















jasc 




END 


















r s 3 


• U 


0AC0CV] Al. 


04C8CV] B 


04D0 CV] 


ORA 




0408 CV] 


D 




C V D 


POI 


0AE8CV] t 


0000 C S J • R 


04 F 8 CV] 


DST 




0500 CV] 


S 




?* CVl 


BP 


?5lCCV] G 


0528CVI A 


0530 CV] 


c. 




0538 CV] 


F 




*0 Cl3 


10 


0000 CS] «I 


K53CCV3 B0T1 


0176 cl: 


11 




0540 CV] 


TOPI 




-.?cl: 


22 


022CCLI S0 


0554 [v: T0P2 


0234 CL] 


31 




055c cv: 


BO T 2 




u cl: 


40 


<?2Ff CLI 50 


0312 CL- 1 60 


056C CV] 


B0T3 


034 » cl: 


61 




7/cv: 


TOPS 


038? cl: 70 


0 4 0? 


cl: 80 


0574 CV] 


TOP 4 


0 4 0 A r l ] 


81 




7? cv: 


bota 


0472 cl: 90 


04rccl: log 


0000 [S] 


• v 










XF.Q 




L 


















Dr 






















^OGRAM LABELS 


« 




















*MA I r, * 


2C0C • V 


2 A 3C 


. COMP 


2C1A PI 




25EC 


. R 


2 B 0 4 


• ZC 


? A 6 c 


*6 


26?C • A 


2 E 6 0 


• Mb's 


29DE »W 




27E6 


EXP 


296C 


A P 


36DE 


alog 


2 A A A • RARG 


2800 


• 5 


2 7E6 AbXP 




2B08 


. ERCNT 


2AC6 


*8 


2908 


.0 


2RDE • U 


3036 
















MTRY-! 


points : 




















?5EC 


.R 


26?C • A 


26DE 


ALOG 


27E6 AEXP 




27E6 


EXP 


296C 


A ^ 


?9D r 


• w 


2 A 3C •CC'NP 


2 A6? 




2 A A4 • R A R G 




2AD6 


$8 


2B00 


. 5 


£P0A 


•ZERO 


2B08 • ERC\T 


2B0 A 


.0 


2B60 • MFS 




2BE A 


.u 


2C10 


. V 



cClA 6 1 



08M0*' -BLOCKS : 
0\E 



NDEPINED SuBROliT INCS : 
O'E 



RtNSEER ADDRESS 20 7 f 

ecuticn begins: 



TCP ASSUMED LARGER THAN BOTTOM 



T r p j = 


e 


2423 


BCT 1 = 


1 


0000 


T C P 1 ° 


p 


2 868 


BO T 1 » 


1 


0500 


T~Pl = 


fr* 


33 78 


BCT 1 = 


1 


000 


T r; P i = 


p 


3960 


BO T 1 = 


1 


1500 


Tqp i = 


0 


<<626 


50^1 = 


A 


2000 


T O p 1 = 


p 


5385 


SOT 1 = 


1 


2500 


Top i = 


? 


6253 


BC T 1*> 


1 


3000 


1! 

CL 

O 


P 


7245 


80 T 1 = 


1 


3 '<00 


Tnpja 


?. 


2 38 2 


BO ” 1 = 


1 

X 


•0 0 0 


T 0°1 = 


? 


9688 


9C T 1 = 


A 


4 500 


T r;Pi = 


1 


1191 


3Q ” 1 = 


1 


5 0 0 0 


T Q P 7 = 


1 


2929 


BC” 1 = 


1 


5500 


T c p i = 


1 


4947 


B? T 1 = 


l 


6 0 00 


T rP - = 


1 


7 3 0 5 


B C T 1 = 




6500 


T c p : = 


2 


<7078 


B0 T 1 = 




7 000 


Tr'P i = 


2 


3366 


3 C ” 1 = 


1 


7500 


TOP 1 a 


2 


3 0 5 


BC” 1 = 


A 


0 0 O 


T C p i = 


3 


8081 


BC” 1 = 


1 


8500 


TqP i = 


3 


7959 


BC ” 1 = 


1 


9000 


II 

<«-{ 

a 

o 

K 


4 


5327 


BC” 1 = 


A 


9500 


T c P i * 


5 


4783 


B0 T 1 = 


? 


0800 


Tr p ps 


0 


0280 


BC”? = 


-0 


0182 


r CP£ = 


3 


04 0 0 


B C ” ? = 




0375 


T C p P = 


2 


^600 


B0”2 = 




0579 


T rpp s 


P 


8 0 0 


B0”? = 




-7 799 


T '' D 2 = 


0 


000 


BC”2 = 




1038 


- 0 p 2 = 




20 0 


B0”P = 


-0 


1301 


T ^ P £> r 


P 


4 7 , 0 


B0”2 = 


- 5* 


1595 


Tnpp = 


7 


16 0 0 


B0”2 = 


- 0 


1 °3? 


T o p ? = 


p 


1800 


BC”2 = 


"0 


2328 


T ^Pp = 


<7 


0 0 0 


B0 T 2 = 


-P 


2 808 


T n P p = • 


p 


220 0 


BC T 2 = 


-P 


t4 1 7 


T fP? = 


p 


2 4 0 0 


30”? = 




4234 


T CP? C 


p 


8600 


BC T 2 = 


-(7 


54 21 


T C p P = 


i? 


2800 


B0”2 = 




7376 


T c p 2= 


e 


3000 


BC”? = 


-I 


1361 


Tc D ? = 


p 


3200 


BC”? = 


-3 


4 815 


t C p 2 = 


0 


3400 


BC”2 = 


1(7, 


62 


T C P 2 = 


0 


3600 


BC”? = 


1 


60 81 


TqP 2 = 


0 


3800 


BC”? = 


P 


8255 


T C P 2 = 


0 


4 0 0 0 


BQ”2 = 


0 


5282 


t 0 p 2 = 


0 


4200 


BC”? = 


0 


3680 


t C p 2 = 


0 


4 4 00 


BC”2 = 


(7 


2655 


T c p ?= 


0 


4600 


BC ”2 = 


0 


1926 


T C D 2 = 


0 


4 800 


B0 T 2 = 


W 


3 7 1 


T 0 p 2 = 


0 


5000 


B 0 ” ? = 


Z 


0925 


TC P 2 = 


? 


520 0 


9C”2 = 


0 


0553 


t CP2 = 


? 


5400 


9C”? = 


0 


0233 


^ c P p e 


0 


5600 


B C ” 2 = 


-0 


0050 


TpPpc 


0 


5800 


90 T 2 = 


-0 


0 303 


t CP2 = 


B 


*00 0 


B0”2 = 




0534 


T C p 2 = 


0 


6200 


BC T 2 = 


7 


0748 


T C P 2 = 


? 


6400 


BC T 2 = 


-z 


7948 


TqPjc 


0 


6600 


BC”2 = 




1136 


t C p 2 = 


0 


6800 


B0 T ? = 


-0 


1315 



D 



T rp 2 = 


0 . 7 0 0 0 


90~2 = 


"0 


•1487 


T0P? = 


0 . 7200 


B0 T ?■= 




• 1651 


t qP ? = 


0 . 7 A ? 0 


90 T 2 = 




• A 9 1 0 


T r P 2 = 


0.7600 


80 T 2 = 


•*P 


• i n64 


T p p p a 


0.7800 


B0 T 2 = 


-2! 


*0114 


T r P ? = 


0 • 8000 


90"2 = 




• ??60 


T C n ? = 


0.8200 


9C~2 = 


-0 


•2403 


T CP? = 


0 . 8400 


9072 = 


-e 


•2643 


t 0 p 2° 


0.8600 


90 T 2 = 




•2681 


T ^pps 


0 .880 0 


90 T ? = 


-<? 


•7*16 


T f"Pp = 


0.0000 


90" r 2 = 


-e 


• 2950 


T QPc' = 


.0800 


90 T 2 = 


-e 


• v 8 1 


T oDp = 


0 . n 4 '• 0 


BC' r 2 = 


w» ^7 


•2211 


VPp = 


0 . Q 600 


9C’2 = 




•3240 


T^Ppc 


0 . c 800 


90’2 = 


- 0 


•3467 


Trppa 


1 . 9000 


90^2 = 




•359^ 



55 



Note that in this case by simple inspection of the 
coordinates of the points to be ploted we can see that the 
intersection will occur outside the acceptable range. 

Note also that in order to obtain the coordinates 
of the points within a range of interest where the intersection 
may occur, different ranges and increments are considered in 
each program. This could be improved if at that time these 
printouts were available and a better idea of the thicknesses 
that should be expected existed. 
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BOTTOM 


assumed 


larger Than 


TOP 




T 0P3 = 


-O . 1776 


BOT 3 = 


l • j- o 0 0 




Top 3 = 


-0 . 1 099 


B O T 3 = 


1 • 0500 




T OP 3 = 


» 0 . 7 310 


R0T3 = 


1 * 1 '/ 0 0 






T 0 P 3 = 


0 . 0596 


B0 T 3 = 


1 • 1 5 0 0 






T C P 3 = 


P. 1635 


B0T3* 


1*2000 






T 0 P 3 = 


0.2782 


30T3 = 


1 *2500 






T 0P3“ 


r . 4GT/2 


BOT 3* 


1 • o -•' 0 0 






T 0 D 3 = 


0 • P 5 0 1 


3973 = 


1 *3570 






T 0'^3 = 


P . 7 0 7 4 


507 3 = 


1 • 4 ? 7 P> 






7 -■ 3 ^ 


^ ♦ 8797 


B 0 T 3 = 


1 • t, S 0 0 






T 0 " 3 “ 


1 « '/fc74 


a o 7 3 = 


1*5-/ 0 0 






' O' 3- 


1 .?71 1 


? 0 " 3 = 


1 * c-^ 5 ” 0 






^ O r- 


1 « <*91 4 


BO'; 3» 


1. % 6 0 /' 0 






y r J -• 

» 


1 *7287 


2- *i *” *Z _ 


‘ ’6577 






~ ~ ~ - - n 


1 • 9836 


r '*. v o -- 

D — 


: > 70 / o 






” 3 j ~ 


2 • 2567 




: * 7" 7 ' 






~ ' j ~ 


2 f 5 A 8 8 


l 7 ~ 3 - 


: • 80 7 -' 7 




r ^ r 3 


? 8 5 9 A 




: • S K 27 




7 * - i: 


3 t l 9v 1 


“ 1 v 7 - 


- , CfS-S / 




~ - 3:,^ 


3 . 5 A 1« 


- •, ' r "j _- 


. >9-' : ‘0 




*' ~ - ^ ■ ■ 


3-9128 


r j ~ 3 - 


? ••/ 0 7 •' 




4 •- 


- v 2 7 9 


r- 3 *' 4 « 


a . 1 4 s 






4 


v ' ■/, 4 0 0 


4 = 


1 > 1 v 7 8 






4 -- 


/ .. 7 6'/ 0 


3 4-4 


• 7 7 






7 ' 4 r - 


0. .•Svr.t 


: ;^ = 


: ’ -333 






7 l =* 


/ , 17-? 2 


7 3 ^ 4 ^ 


0.0061 






- 4 ~ 


/ . 12- .. 


3 3 * =» 


0 j 9f89 






4 * 


'• .14/,' 


£ ; * r 4 - 


0 '9217 








- * 160-' 


5 ; ‘ 4 = 


0-8845 






7 ' 4, ~ 


' . i B'.'O 


£; 4 -= 


0 - 8 4 7 3 






4 - 


* - ? •-' 


3 3 4 = 


0 8 1 M 1 






7 - 4 -■ 


.22 '• 


L ^ ‘ 4 ii 


• 7 7 2 8 






T ~ r c ' 


-24 


•• ; " v - 


3 7 35 4 






7 H • 


n ? /. * 


r wj. ~ 


M 69 84 






■ 


- P ^ 1 


* 3 ‘ = 


0*6612 






T* 30 ^ ~ 


' 3 





7 6 2 






~ 2 7 V = 




' ~ c = 


0 - u k 6 8 






~ ' - .. = 


7 * 3 *+ ’7 


r* *> ” i = 


0*6496 






3~<- 


. -363-? 


5 0 " «• = 


0 - 5 1 2 4 








„ o 


33-- = 


0-4752 






*"7-4- 


- 'k*3 . 


3 - - a., s: 


0-4379 






^ ~ c — 


^ - 4 2 T ' 


03 “> *“ tx — . 


0-4707 






~ “* - z = 


/ - 4 >+ r ~ 


: 3 T ^ = 


0'3.35 






•“ ■ ^ ... 


* U •- 


P ,1 * 4. = 


0-0263 




- «.v = 


3 - 4 : 3 J 


3 7 — = 


! V 0 ? Q | 




™ •> O' ^ _ 


- . rr ^ 

•-- V 


5 3~u.- 


0-2519 




-u. ” 


2-523' 




Z ' p 1 A 7 




" ^ = 


* -54-3 


33 ^ " 


0-1775 




" ^ _ 


0 . 5 4 7 2 


P'“'~/L8 — 


0 - 14 03 




3 *-■ ^- = 


^ r* ^ o 


5C”4 = 


7-1330 




• ^ n _ 


/ * f 7 7 7 


p — — 


? a 71 < 5 ? 




3 7 = 




2 1 v ^ — 


3 - - 3 ? 6 




* 7 - /* = 


7 * - 4* 7 7 


33“^ = 


M 7 • ^ p A 




* 3 - 4 = 


2 * : o :; 7 


3 3^4 = 


•• 3 * ^ u ^5 ? 




"7 -4.= 


r • 6 > r 7 


3 3 " U rz 


- 6 -3-3 Z 1 
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*”* o c • r~ 


r* 


. 7 " ;• 7 


S C = 


» 7 « < T 7 P 




7 


*7^7? 


7 7 ~ ^ = 


• 7 • 17 7^ 


~ c ° = 


r> 


• 7 ^ 7 7 




*■ 7 • t ^ £ 




fr 


* 7 A 7 7 


^'/Ar 


-7 • ?■? 1 9 


~ c r - 


</ 


.yew 


3 ^ ~ /- r: 


-7 » ?<: 91 


~ ^ lL = 


r> 


m < r>~? 


p 7 ~ *. = 


-7 - 3"63 


" - ^ 4. = 


o 


* c P 7 '/! 




- 7 • ? t R 5 


" Q ^ — 


7 


. - /j 7 7 


r\ — 


-» £ * ^ C Cf 7 


'C 3 *- 


7 


♦ r /Ss 7 7 


o o — = 


- 7 * v 1 7 Q 


~ r - /.• = 


rt 


. v c 7 7 


p c * = 


-7 • 4.^51 


~ ^ 5 * = 


7 


. r iff r> ,£ 


P 7 “ U rr 


ap? 


m :^,= 


r 


• o *? ~ r> 


^ r> - u 


<*• c «c:';Of 


~ r - l = 


/y 


* c 4 - '? 


^ ^ " /-•- s 


-o r.P 


"^/i = 


rj 


. -.A 7 7 


: ^ * L - 


— 7 • 4 f * <4 7 


* r ^ 4 = 


f? 


.097/ 


co- /-. — 


-? • i4 12 




1 


. r 0* '? ;> 


^ C ” y ~ = 


- 7 • A 7 ? 4 



Cf vio OPERATING SYSTEM VERSION ?. REVISION 01? 03/04/74 GENERATED 
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/ 


jnp 
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LINES 


PRINTED PR? 
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/ 


K LN 
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67 
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READ 
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/ 


K cn 
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00 
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/ 
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17 


M ODEL 


70 SECONDS 


*25.00 


/ 
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00 


MODEL 


50 SECONDS 


*1 ?,53 


/ 


u 0 U R 


• ?7 








total 


CHARGE * 


•87 



:RFEIR 490 14731 LOGGED OUT 75/06/74 15:25- 



* 12.42 LEFT AFTER 



05/06/74 



19 LOGINS 
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SUMMARY 



Length 


— 


o 

o 


Beam 


= 


AL/8 


Draft 


= 


B/3 


Depth 


= 


AL/9 


DST 


= 


20000 psi 


S 


= 


120" 


Material 


- 


STEEL 


TOP 


= 


0 . 03 " 


BOT 


= 


1.14" 


STH 


= 


19771.88 


SBH 


= 


9579.87 


SBCRIT 


= 


9788.37 


DITOP 




37.42' 


DIBOT 


= 


18.13' 



Note in this case the large value chosen for S and 
associated with it the large value required for the thickness 
at the bottom to avoid buckling which nevertheless does not re- 
sist stresses larger than 9788 psi in compression. 



62 



If we look now at these three cases, all for the same 
material and. set of dimensions (same bending moment) we can 
see how the thicknesses and neutral axis position changed as 
the value of S was changed. 



CASE 


s 

inches 


TOP 

inches 


B0T 

inches 


SBH 

psi 


SBCRIT 

psi 


DIT0P 

ft 


DIB0T 

ft 


1 


10 


0.468 


0 . 166 


30000 


30000 


22.06 


33.05 


2 


30 


0.4125 


0.4125 


20000 


20000 


27.78 


27.78 


3 


120 


0.03 


1.14 


9700 


9700 


37.43 


18. 13 



In this table, we can more readily make conclusions on 
the influence of the value chosen for S on the neutral axis 
location, and also a quick look at the resulting values for 
the thicknesses may give an idea of the weight associated with 
each case. 

In the following cases the dimension parameters are 
changed (increased) and consequently the bending moment and the 
thicknesses of plating will be higher. 
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STRENGWT OF SHIPS-QN ThE NEUTRAl AXIS LOCATION 
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T O c 1 = 


0.0844 


T np 1 o 


0*1554 


Top i = 


0 • ? 6 7 3 


T 0 P 1* 


0.4414 


T c < = 


0.7145 


T0 D 1 = 


1 • 1570 


T 0 P i = 


1.9240 


T 1 p < — 


3.4364 


T 0 ° 1 = 


7*3849 


II 

^ c 

a 

O 

K 


36 .2430 


j. 


-20.76*8 


T n D <i — 


-9.8191 


” c p i = 


-7.13*7 


”0P 1 B 


-5*9771 


T o c : = 


-5*3601 


T P o , - 


-5.000 0 


TqBi z 


-4.7*16 


t op i = 


-4 . 14498 


T ^ p l = 


- 4 . r; 7 4 7 


ToD< = 


“4.63?* 


^ r>p 4 = 


-4.4313 


T o p i = 


—4*5450 


T r>D^ — 


-4,5748 


T cj D 1 = 


-4.6173 


T oPl = 


-4.4700 


T C P 1 = 


-4 . 7309 


T n D < - 


-4 . 79*6 


T C P 1 = 


-4 . *,720 


T n P i s 


-4.9502 


T~Pl *= 


-5* 03?5 


t^d-_ 


- 5 . 1 1 s 3 


T 0 p 1 = 


-5.2071 


T np < s 


-5.29*6 


T O P 1 b 


-5 • 3924 


TnP <: 


-5. 4**3 


T C p l a 


-5.5860 


T 0 p 2 = 


1 .0000 


T 1P 2 = 


1 .5000 


T }P£ = 


2.0002 


T ^ P J = 


2 . 5000 


T ^ D p c 


3.0000 


T 0 P 2 = 


3.5000 


T 3 P 2 = 


4 .0000 


t 0 p 2“ 


4.5000 


t 0 p 2 = 


5 • 0000 


T nPp = 


5.5002 


”0 D 2 = 


4 .0000 


T 0 Pp = 


6*5000 


T ^Dpa 


7.0000 


T OPpB 


7.5000 


T 0 P 2' S 


*.0000 


T o p a* 


8 . 5 7 0 0 



BQT 1 = 


0* 1000 


B0T1 = 


0* 1500 


B0T1- 


P«?000 


3QT 1 3 


0 • 2500 


bo” 1 = 


0 .9000 


B0”1 = 


0 * 3500 


90” 1 = 


0*4000 


90” 1 “ 


0 *4500 


90” i = 


0*5000 


BQ T 1 = 


0*5500 


BO” 1 = 


0 *4000 


90”1 = 


0 • 6500 


BO” 1 = 


0*7000 


B 0 T 1 = 


0 * 7400 


11 

T-t 

* 

O 

cr 


0*8000 


4 n - 1 _ 


0 • 8 5 0 0 


BO”?. = 


0 • 9 0 0 


bo t i = 


0 . qc00 


9 q - 1 - 


1*0000 


BO” 1 = 


1 *0400 


30” 1 = 


1 • 1 070 


9C ” 1 = 


1 * 1 500 


B C ” ?. = 


1 • 2000 


90*1 = 


1 * p 4 0 0 


B 0 T 1 = 


1 * 3000 


BC ” 1 = 


1 * 3500 


90” 1 - 


1*4000 


30” 1 * 


1 *4500 


90” 1 » 


1 *5000 


BO” 1 = 


1 *54^0 


90- 1 s 


1*6000 


90- 1 z: 


1*4500 


90” 1 = 


1*7000 


90” 1 s 


1*7500 


BO” 1 = 


1*8000 


BC T 1 = 


1 • 8500 


90” 1 = 


1*9700 


90” 1 = 


1 • 9500 


90” 1 = 


2*0000 


B0’2 = 


-0*4473 


B0”2 = 


-0*7727 


B0”2 = 


- 1 * 1 *29 


B0”2 = 


“1*8315 


B0”? = 


-3 • ?46 1 


90”2 = 


“11 *8?71 


BO” 2 = 


5*8*10 


B0”2* 


1*8850 


Be ” 2 = 


0 * * 1 27 


90”?3 


0*2578 


Bo”2 = 


-0*1140 


B0”2 = 


-0 • 4000 


B0”2 = 


-0 *63*9 


B0”2 = 


-0 • 8492 


BC ”2 = 


- 1 *0408 


B0”2 = 


-1*2195 



Tf'Ppa 


9 . 0 0 0 0 


T nPpe 


9 . 5000 


T C?P 2 c 


1 0 . / 0 /, 0 


T r P £ c 


1 0 . P 0 0 0 


T(>Ppsa 


1 1 . 0000 


Tnp ?c 


1 1 .^000 


T c p ?= 


12.00/0 


T p £ — 


12 * bP'/M 


T o c v = 


13 * 000 0 


top^b 


1 3 . b 0 0 0 


TfiPoa 


14 . 0 0.0 0 


t Q p ? = 


14.50./0 


T r P p = 


15*0002 



B0 T 2 = -1*3*90 

BQ t 2 b - 1 • 55 1 6 
BC T 2 = - 1 - 7.790 

9 0 T 2 b - 1 • 5 ? 6 ? 3 
BO T ?= -2*0124 
B0 T 2 = -2*1599 

BC T 2“ -2*3/52 

RO T ? = “2*4^27 

BO T P= -2*5907 
B0 T 2= -2*7914 

BCT2- ”2 • K 7 1 1 
B0 T 2 = -3*0/98 

5C T 2= -3*1477 



6? 




bottom thikness (in.) 



68 



OnTTD 


M ASSUMf'O 


LARGER THAN 


T 0 p 


T n P ? » 


"• 0 • ■" 3 1 6 


RO T 3 = 


(? • 1^00 


T ^ P ;j s 


-0 . °2p7 


30 T 3 = 


(1 . ^ 5 0 0 
('•9000 


TriPo* 


- C • - 0 6 6 


RO T 3 = 


T <*> P 3 = 


4 . 4 4 9 


?0 T 3 = 


0*2400 


T ~ D 3 = 


0 ♦ 1 3 1 8 


^0 ~ 3 = 


0 . 3**00 


T - p 3 = 


0 * ? 6 1 3 


3C~3 = 


0*3500 


’'Oje 


0 t 4 4 v* 2 


BO ~ 3 = 


e • 4 0 0 0 


~'-P? = 


?■ . 67b 8 


R 0 “ 3 = 


0 • A K 0 0 


w 3 = 


••” * 9 7 - 2 


3q - 3 - 


0*5/00 


~ r . ° 3 = 


1-3455 


=• r " 3 = 


j*5~00 


“ ^ 3 = 


1 - 7906 


BC“3 = 


0 - 6 ? 0 2 


* ' ^ ^ s 


3 - 3?6 ? 


3 C ~ 3 = 


0*6 5 ? ? 


T ~ ^ 2 = 


? - 9b?'/ 


w r - 3 - 


7*7000 


T pcs 


3 - 6 7 6 4 


D,rrP- 


7-7500 


T '* >' ? - 


4 - * (?. ; 1 


O r« ~ ~ _ 
• — — 


7 - p / 0 9 


"rP *.- = 


b' - A 5 1 1 


3C"3.t 


7 - P50Z 


c 2 = 


6 - 5 1 5 5 


2 . r* — c .. 


*> »* O * C* 2 


" p 


7 - 7 1 * 7 4 


EC~3 = 


/ * 9^72 




9'7?o9 


P 7 ~ 3 5- 


1 * 0/00 


“ ~ ^ 2 sr 


1 '/ - r -021 


*. — ~ . 


1 * 0-00 


T ^ r- 2 = 


13-11 ;C 


3 C “ 3 = 


1 * 1 0 00 




1 3 ■ ;- 9 t 7 


E C ’ 3 = 


1 * 1500 


7 3 - 3 = 


lb- op 73 


^;p = 


1 * 2 0 0 0 




1 7 • 333'? 


2 , 0 — c _ 


1 - 2-0 0 




? '7- PI 57 


^ ^ - -7 _ 


1 - 0 *> 7 , 0 


T p 3 « 


?? - '•??<!- 


C* r' ** ^ __ 
- _/ *— — 


1 * 0-00 




25 - 34/7 


? 0 “ 3 = 


1*4/00 


7 7 ~ 3 = 


? ?' - 197? 


7 C~? = 


1*4370 


~ ~ p c = 

> 


31 -”931 




1 *5**00 


^ y ^ — 


34 . c 3 3 4 


EC'3 = 


1 * 4 *- 0 0 - 


T r P *: r 


3? - P74 


O, ^ ~ _ 


* - r-'-ZZ 


“ ^ ^ p = 


4 1-74 .-; C* 


*2 r* 1 _ 


1*4900 


*^-7“ 


4 <= - 7 7 p 4 


D^^p- 


1*700 


” ^ P ~ r 


4 9 - ?°76 


C, H> ” *2 = 


1 * 74^0 


^ P X 


[3 4 . 0 4 *7 7 


2 ; r* — *3 _ 


1 - p '*00 


* « n * 

r* - — 


co . 3 f c: 


O, p - O - 


- • p r -; 7 7 


T ^PP = 


6 0 . 0 a a ot 


Rn^3 = 


1 • 0 7 72 


T n°ps 


9 9 . P P P 3 


PO T Ps 




T nDp = 


74 . 73-09 


30~? = 


2 . 


"nP 3 = 


R 0 . c ? p 8 


R0"3 = 


2 • ^^07) 


^ 0 D pa 


S4 . c 0?2 


BC~3 = 


2 • 1 00^5 


T pD ? = 


9? . 979? 


R2’ r 3si 


2* *500 


" r C D 3 = 


99*6604 


B0^3 = 


2 • ?000 


* C D 3 = 


1 PA . ^534 


R 0 ~?« 


2*2500 


-QP ? = 


1 1 3. 96 52 


30^3 = 


2*3000 


t ^PT = 


1 ? 1 .6029 


RC^3 = 


2 * 3 b ?0 


T -'0-3 = 


1 ?9 *9735 


30^3 = 


2*4000 


T -ois 


1 3 7 • 0 8 s 0 


BC^3 = 


?• 4500 


TqP 3 = 


146.5416 


90 t 3* 


2 • c 0 0 0 


" r C D 3 = 


195 • 5594 


30^3 = 


2*5500 


O D 2 | t= 


164 . 92f 5 


R 0 ^ 3 c 


? • 6700 


t : c 3 c 


174.4679 


BC^3 = 


2*6500 


t ''PTs 


1 R 4 . 7 p 4 b 


9C’3“ 


2 • 7000 


7 Q w 3 c 


1 95 . 2838 


80^3 = 


2 • 7500 


’C p 3 = 


206 • 1 7p6 


90’3* 


2*8000 
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T pP3 = 


217 .45*1 


B0 T 3 = 


2 


•8500 


T p P 3 c 


229. 1475 


50 t 3 = 


? 


♦ 900 0 


T c P 3 = 


24 1 .2473 


B0 t 3* 


? 


• 9500 


^ lP 3 n 


253.7654 


BC t 3 = 


3 


*0000 


r o p 4 = 


1 . :* 00 M 


B0 T 4 = 


16 


• 1 639 


T n P ^ s 


1 .50/0 


BO T 4 = 


14 


• 2342 


V^/iC 


2 . 0000 


B0 T 4 = 


1 ? 


.3044 


T n o u B 


2.5000 


B0 T 4 = 


10 


•3746 




3 • 0000 


B0 - '4 = 


8 


• 4449 


T O p 4 = 


3.^000 


B0 T 4 s 


6 


•5150 


T ^/ + = 


4 . 0000 


B0"4 = 


4 


•5*52 


T ^P4 = 


4.5000 


B0 T 4 = 


? 


• f. 554 


T ' P /* = 


5 • 000/ 


BO ” 4 C 


0 


•7256 


T nP4 = 


5 . 5 0 0 0 


00 T 4 = 


-1 


• ? 0 4 2 


T ^P(i = 


6 • / 0 v Z 


BO - ' 4 s 


-3 


• 1 340 


T p P 4 = 


6.50O0 


B0 - '4 = 


""5 


•0a38 


Tf 'P4c 


7 .'*000 


B0 T 4 = 


- 6 


• 9936 


7 r* D 4 — 


7.5000 


B0 T 4t= 


~ 8 


• 9235 


T ^ P 4 = 


R • 0 000 


gn*4 = 


-10 


• «532 


T -* P 4 = 


8 . 500 0 


Rn' r 4 = 


-IP 


•7829 




9.0000 


bo t/ * = 


-14 


•7127 


11 

-J 

a 


9 . k 0 0 0 


BC~ 4 = 


- 1 6 


• 6A?b 


T n P 4 = 


10 .'•*000 


BC"^ = 




•5723 


T p u = 


10.5000 


gn T 4- 


-20 


• 5 0 2 1 


T i' P a = 


1 1 . 0000 


B0 T4c 


-2? 


•4319 


T c 3 - = 


1 1 • 5000 


pn T 4 = 


-24 


•3617 


T e p A = 


1 2 . 0000 


B0' r4 = 


-26 


• 2° 1 5 


t = 


12.5000 


RO T 4 = 


-28 


• 5213 


T ;P<* = 


13 . 0000 


B0"4e= 


-30 


• 1511 


T ~D^ = ' 


13*5000 


BC^ = 


-3? 


• 0908 


T 


14.0000 


3n T 4 = 


— 3 4 


• 0 1 0 6 


TlD t : 


14.5000 


BC ^ 4 = 


-35 


• 9 4 0 4 


T 0 p A = 


15.0000 


90^4 = 


-37 


• 8702 



Vi: OPERATING SYSTEM .'ErSION 1 REVISION 012 03/04/74 GENERATED 



■Ob 


U ANDL I NG 


CHa 3 GE 


* .35 


/ 


JOB 


53 


LINES 


PR 


InTEO p R2 


* 1.25 


/ 


K LN 


67 


CAROS 


PEA- 


s 1.5/ 


/ 


K CD 


0? 


c lotte 


R 


VICTORS 


s .25 


/ 


1/07 


17 


v ODEL 


70 


cECONDS 


s?5.0/ 


/ 


MO'JR 


00 


w odel 


9? 


SECONDS 


M2. 5/ 


/ 


UOIJR 










TOTAL 


C-iARGE 



UP 14731 LODGED '*,UT 7 5/06/74 



•35 
•3? 
. 10 
• 00 
• 12 
• 00 

* • 89 



17126. $ 9.62 LEFT AFTER 



05/ 06/ 7 A 0 



23 LOGINS. 
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p-pa 

c pa 
e.TC 

s l? 

t’ris 
Z '2 k 
Z '3«' 
Z 3« ? 
P 3*’ 
Z 3? 
a ' 3 8 
Z -38 

a *c 

< -c 

a- b* 

z 

a 'sC 
Z pC 

a *8 
a r A 
a 5 *» 
a 5 a 
Z 3C 
a'cA 
a’rp 
a- ac 
a* *c 
a s 5? 
a- 7 a 
a:sE 
a b a 
a- oi 



c *AL* IS the lenght 
AL* 10Oe. 

B-AL/5.75 

c * OK A » IS THE ORAFT 

D k A=B/3* 3 
O s AL/1 A • 

c *POI* IS THE POISSON'S RATIO OF THr MATERIAL 

c for STEEL 

p 0 I = 0 • 3 

c * E » is The YOUNG'S MO n ULuS OF tHE MATERIAL 

C FOR STEEL 

E e 30»-10«*»6 

C .OST* IS THE OESIGn S t REsS 

C)S 1*18000. 

•S. IS The SPACING OF LONGITUDINAL STIFFENING MEMBERS 
S = 30 • 

c .Bn. is the Rending m 0 ment 

PMb(U ,i DRA*aL*aL*0 , 7s) /(3S..35. ) 

G c e. 30396 *S*S* ( 1 . -POI *POl ) /£ 
c .TOP- IS THE THlCKMESq OF the top PLATING 

c -POT- IS THE THICKNESS OF THE BOTTOM PLATING 

jOpcS • 1 P 
f3OT*0*62 

A c 5 60 • *Bm* (d*POT + 0* TOp ) 

[)En = 0* ( 3. *B*B*TOP*ROT + E. .B*0*T0P* ( TOP + BOT ) ♦ ( 0*0* TOP* TOP ) ) 
C * S TH * IS ThE ST r ESS On ToP OUE to HOGGING 

STH = A/[)tN 
W«ITE(S/60 )STh 

60 FORMAT! ' ' ,'dV\y, 'STHa'plH.3) 

C=560. -f'M-(B*TOP + O*T0p) 

C • SRh * IS THE STRESS On b 0T DUE T 0 HOGGING 

sBh=c/den 



a:C6 

z: 12 
a : pc 

I'.rC 

i?ZZ 

i’ 2 ? 



W«ITE(5»90 )SB W 

90 FORMAT!' ' , 2/iy, ' SBH= ' r 1 0 . 3 ) 

C * SBCR I T " IS TnE CRITICAL STRESS ON BOT 

SBCRIT=BOT-BOT/G 
WHITE (5# lit*. )SPCPIT 
lit? FORMAT!' ',17y, 'SBCRl T = ' F 10 . 3 ) 



a ?62 
a?6S 
a?« c 
a?** 
azc*. 



p*E.*D-Tof > + B*(TOP + POT ) 

C * 0 1 T 0 P » IS THF DlSTAN r E OF THE TOP FROM THE NEUTRAL AXIS 

DIT0P = 0*(B»B0t+G-T0P) / P 
R«ITE(5/ A^jOItoP 

Ap FORMAT)' ' , I8x, 'OlToPs'FlO. A ) 

C .DIBOT- IS THF DlSTANrE OF THE BOTTOM FROM THE NEUTRAL Ayis 



i 5 CA DlBOT=D“ ( 6*T0o+0«T0P ) /P 

wRlTE(5>5(5)0ln0T 

5c* FORMAT!’ S 18x> 'DlBOT='FlO.A ) 
2 ^ 2 1 £ N 0 



aa csi 


•U 


032a CV] 


AL 


03?C C V I 


B 


033 A CV] 


ORA 


033C CV] 


D 


♦ a CV3 


Pol 


e3Ac cv] 


E 


000? Cp] 


. R 


03S C CV] 


OST 


036 A CV- 


S 


i- CV] 




037C CV] 


G 


0388 Cv3 


TOP 


03^0 CV] 


BOT 


0398 CV] 


A 


* s CV] 


Den 


03BC CV] 


STH 


01 7 a C L ] 


60 


0000 CS] 


>9 1 


03C0 CV] 


C 
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•KV] SqH « 1 E 2 CL 3 90 03CRCv3 SBCRIT F 2 2 8 C L J 110 03CC CV) r 

^(V3 0 1 TOp ?2A8CU 40 03D4C V 3 DlPOT 03F4CL3 50 0000CS] .V 

*EO L 

)F 



tOGRAM laPELs • 



? <’ 7 0 


• u A In* 


2A68 


. V 


PR98 • ( -OMp 


2 A 7 6 


pi 


24 48 


.R 


2960 


• 2 


2 5BE 


*6 


2488 


. A 


29BC «yES 


2°3A 


. W 


264? 


EXP 


27C8 


A I 


2534 

2964 


A log 

•0 


2900 
2 A 3 A 


• R A R G 

• U 


?95C • 5 
3B92 


2642 


A EXP 


2964 


. ercnt 


2932 





’RY-POlNTs ! 



2448 


• P 


2488 .A 


253 A fl, OG 


26^2 


AEXP 


264? 


EXP 


27C8 


A I 


?83A 


• W 


2898 .COMP 


2«PF 5 6 


2900 


.RARG 


293? 


*8 


29bC 


• 5 


2960 


•?EPo 


2964 . ERCNT 


2966 • 0 


29t3C 


. MES 


2 A 4 6 


• U 


2A6C 


• \} 


2 A 76 


I 



















'•mon-Blocks ; 

\ r 



.EFINEo SuBPOu 1 INES : 
'•E 



ANSFER ADDRESS 2070 

cut ion bfginS: 

STH= 18085,023 
SBH= 47973.711 
SBCRIT«= 4632^.625 
D I T OP= 19.5551 
0 I BO T = 51.8734 

End 



F Vio 0 °FPAtInG SYSTEM version 1 R r VISION 012 03/04/74 GENERATED 05/06/74 



JOB 


HANDLING CHARGE 


* • 35 


/ 


JOB 


* 35 


85 


lkes printed pri 


t 1 ,?5 


/ 


K Lv[ 


• 1 1 


5** 


ca=ds reao 


* 1*^0 


/ 


K Cq 


• 0 K 


fC 


PLOTTER VECTORS 


* • ?5 


/ 


1000 


• 0 ^ 


13 


KOdfL 7,7 SECONDS 


*25 . ^0 


/ 


HOUR 


• 09 


0 ? 


,v, 0n EL 87 SECONDS 


512.52 


/ 


‘-04JP 


t 






total 


C H A P 0 f ® 


• 63 



REIR 4qa j 4 7 3 1 LOGGEO OUT 7,5/06/7,. 17!36. * 3«11 LEFT AFTER ?5 LOGINS* 



SUMMARY 



Length 


= 


1000 ft 


Beam 


= 


al / 5.75 


Draft 


= 


B/3.3 


Depth 


= 


AL/l4 


DST 


= 


18000 psi 


S 




30 " 


Material 


- 


STEEL 


TOP 


= 


5.12" 


BOT 


= 


0 . 62 " 


STH 


= 


18085. 02 


SBH 


= 


47973.71 


SBCRIT 


= 


46323.82 


DITOP 


= 


19.55' 


DIBOT 


= 


51.87' 



We might say here that the value chosen for S is small 
and this accounts for the compression at the bottom to reach 
such a high value. 

In the next case, the value chosen for S is increased to 



60 inches. 
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???* 

iV ?4 

J-’.’C 
J.'lR 
J 18 
J 2 fc 
30 
J 3«’ 
^32 

i*. 38 

.’ 73 ? 

38 

; *c 
i? <»c 

5* 

54 
.VbC 
: 5C 
: ’ ? s 

iVr* 

OC 
i r “ 
- 1 EC 



c 

c 

c 

c 

c 

c 

c 

c 

c 



c 



* A|_ * IS THE LENGHT 

aL= i 000 • 

B = AL/5 . 75 

»ORA* IS THE oRAFT 
D r A*b/3 , 3 
0 K AL/1 4 • 

* PQ I * IS THE POISSON'S PATIO OF THE MATERIAL 
f Ok STEEL 

pOI =0 • 3 

*E« IS THE YOUNG'S M.OqULUS OF THE MATERIAL 

FOR STEEL 

E=3P.*10.**fe 

* OST * IS THE DESIGN StRESS 
qST = 1 8 0 , v 3 • 

*S» IS The SPACING op LONGITUDINAL STIFFENING MEMBERS 
S = 60 • 

* Bn * is the Bonding M 0 MtNT 

BH=(B»DRa«AL*aL* v1, 7= ) /( 35 •*’35. ) 

G=0* 30396 *S*S* < 1 . -PO I *P0 I )/E 

A s 9 • » B * P • 0/ ( 1 /* 0 • * U.m ) 

C = 3 * * B » D * D / ( 7 p . * B M ( 

F = 3**D»0*D/( 1 i , (1 • • ti M ) 

THICKNESS OF T O p ASSU m ED LARGER THaN BOTTOM THICKNESS 



»!2C 

; i ,i 

1 6E 

’ 16 E 

fct 

■»176 
! eE 
1 E 2 



rjPlTE(5/ 12) 

IP FORMAT (' 1 25x, ' TOP ASSUMED LARGER THAN BOTTOM') 
C * T OP * IS THE jHlCKNESq OF THE TOP PLATING 

*00 T * IS THE THICKNESS OF THE BOTTOM PLATING 
B 0 T 1 = 0.02 

ll tOP 1=B»60T 1 *B0T1 *BOT 1 / ( <B+D)*G*OST-D*BOTl*pOTl) 
wRITE( 5 * 20 )T 0 pi>B 0 T 1 

FORMAT)' ' / 25 X / ' TOP* = t F9 . 4, 5X, 'BOT * = ' F9» 4 ) 







t)OTl«=BOT 1 0 • 0 p 


’716 




IF (B0T1 • GT . 2 • ) GO TO 3* 


:=>22 




GO TO 11 


' 5 2 C 


3? 


T0P2=1 . 


* ? 3* 


31 


B0T2 = T0P2* < 12.*0- (C + F , *0ST * T0P2 » / < DST * ( A + C ) * TOP?- 12 


Jr??? 




W«I Tt < 5 # 40 » T0P2, GOT? 


‘254 


40 


FORMAT) ' '/25 X , 'T0Pp*»F9.4,5X, 'B0T2='F9»4 ) 


1 : 2 D C 




T 0P2 = T 0P2 + 0 • 2 






I F ( T 0P2 • GT . 6 . ) GO TO 50 


' " r A 




G° TO 31 


'?r[ 


50 


CONT INuE 




C 


THICKNESS OF BOTTOM ASSUMED LARGER THAN TOP 


! ! J rL 




WRITE) 5,60 ) 




60 


FORMAT (' 1 25X/ ' BOTTOM ASSUMED LARGER THAN TOP') 






BCT3 = 0 . 4 




61 


T0P3= ( (EUD)*B0T3*90T3.BUT3-D«'G*DST-vB0T3)/{p*G*DST) 


nsk 




wRITE(5>70)TOp3*BOT3 




7? 


format ( ' ' , 25x , ' T OP 3 * 1 F 5 » 4/ 5X > » BOTt* ' F9»4 ) 


►2E8 




bOT3 = BOT3 + 0 * 0 ? 


>-£C 




IF ( BOT 3 *GT . * 85 )G0 TO *0 


12 * E 


80 


yOPA * 1 , 
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.-.P6 81 0 OTAb( 12. * ( U*f> )- ( A + c ) *OSt*TOP<* )/( ( C + F ) *OST ) 

.i(,A wRITE(5>90)T0oA/B0Ta 

, fc6 E 90 FORMATf' ',25 x> ' T0Pa=iF9.4,5X, <B0Ta='F9*4 ) 

:i .g6 TOPA*TOP<f»0» 2 

. a2 iF ( TOP<* .GT.S.s )CiO Tn i0« 

! 1 0 <* GO TO R1 

wB? i P10 c 0n TINijE 



■*. b8 
. CS3 


•u 


end 

04BC CV] 


AL 


0 A C 4 C v 3 


B 


0 4 C C C V ] 


DRA 


0404 CV] 


D 


C tv] 


POI 


04E4 CV] 


E 


000.0 CpJ 


.R 


04F 4 CV] 


DST 


0 4 F C C V i 


S 


• CV] 


3 m 


0518 CV] 


G 


052a C v 3 


A 


05?0 CV] 


C 


053C CV] 


F 


cu 


in 


0000 i C S3 


(» I 


C’bA0 C v 3 


B0T1 


0176 CL] 


1 1 


0548 CV] 


TOPI 


? CLI 


2? 


022C CL] 


30 


0558 C \/ ] 


T0P2 


0294 CL] 


31 


055C CV] 


BO T? 


A Cl] 


<* j» 


02FE CL] 


50 


0312 C l j 


60 


0574 CV] 


B0T3 


0348 CL] 


61 


r CV1 


TOP? 


03B8 CL] 


70 


0 3 F E C L ] 


80 


054 4 CV] 


T0P4 


0406 CL] 


81 


• CV} 


BOTA 


0A6E CL] 


90 


04 R p C L ] 


100 


0000 CS] 


.V 







• LO L 



~GR A w LAPELS : 


?,'70 *ma I m* 


2C20 


.V 


2A50 


• fOBp 


?A7 b Sft 


2640 


. A 


2B 7 A 


» mL S 


cfcF2 AlOG 


2AB8 


.PARG 


2B 1 4 


• 5 


‘B1C *0 


2BF 2 


.U 


3D4 A 




’Ry-points : 


2600 » R 


2640 


. A 


26F2 


A 1 OG 


? 9F 2 -w 


2 A 50 


.COBP 


?. \1 b 


* 6 


2818 • 7 EP 0 
5 C2E Pi 


2B1C 


« ercnt 


2B1E 


•0 



2C2F 


f I 


2600 


. R 


PB 1 8 


• 7 


29F? 


. w 


27FA 


EXP 


?9«0 


A T 


27 F A 


A EXP 


2B1C 


•ercnt 


2 AE A 





27FA 


aEXP 


27F A 


EXP 


29«0 


A I 


2 A 58 


.RARG 


2 AE A 


5 8 


2B 1 4 


• K 


2B 74 


, NES 


2BFE 


• U 


2C24 


• V 



*-.ion-bi_ocks : 
*»£ 



■cFineo Subroutines: 
' E 



s.NSFER ACDFESS 2070 
COT I ON BFgInS: 



T 0 P 1 b 


0 


.7/(76 


top i = 


0 


» 70 1 2 


TOPls 


0 


.//pi 


TOP 1 a 


0 


.77q3 


T 0 P 1 = 


0 


.77.49 


TOPIb 


0 


.7777 


TOP 1 = 


0 


• 7 0 q 7 


TOP 1 = 


0 


• 7 1 p9 


TOPI* 


0 


.71*9 


top i = 


0 


.7816 


TCP 1 c 


0 


.7271 


TOPIb 


0 


• 73-55 


TCP 1 = 


0 


. 7 4 n 9 


TOP 1 = 


0 


• 7 . 4 q 4 


T 0 Pi = 


0 


. 7 Sq 1 


TOPls 


0 


.7/7/ 


TOP 1 = 


0 


. 78 p 4 


TOP 1 = 


0 


.09*2 


top i = 


0 


•Ul 6 


TOP 1 = 


0 


• 12*8 


top i = 


0 


. 1^78 


TOP 1 = 


0 


. l 6p9 


T OP 1 n 


0 


• 1 9p 1 


TOPIb 


0 


•2178 


T OP 1 = 


0 


• 24 5 9 


top i = 


0 


•?7*9 


TOP 1 * 


0 


. 3l 78 


T 0 P 1 = 


0 


•34 7 9 


TOP 1 = 


0 


• p 8 q 6 


TOPla 


0 


. 4 3 


TOPIb 


0 


» 4 8 j 6 


TOP 1 « 


0 


• 53 48 


TOPl« 


0 


• 59p5 


TOPI, 


0 


. 6 B 50 


TOP 1 = 


0 


.7248 


TOP 1 a 


0 


• 8 7 7 3 


TOPI, 


0 


• 8 8 p8 


TOP 1 = 


0 


. 9 7 q 0 


T 0 P 1 = 


1 


.7.7^ 


TOPI* 


1 


• 1 7 q 9 


Top i = 


1 


•29p5 


topi. 


1 


• 4 2 p 9 


TOP 1 * 


1 


• 57 ? 4 


TOPI. 


1 


•7376 


TOPI* 


1 


. 00^4 


T OP 1 = 


2 


.79o0 


TOPI, 


2 


• 3 1 4 1 


T 0 P 1 * 


2 


. 5bq7 


TOPI, 


2 


. 8 2 | 8 


TOPI, 


3 


• 1 2?9 


TOP 1 * 


3 


. 46p6 


TOPI, 


3 


• P 4^7 



0.7270 
0.0400 
0.0600 
0.0300 
0 . 1 000 
0.1200 
0.1400 
0 .1600 
0 . 1 800 
0.2000 
0.2200 
0.2400 
0*2600 
0.2800 
0 . 3000 
0.3200 
0.3400 
0.3600 
0.3800 
0.4000 
0 . 4200 
0.4400 
0.4600 
0.4800 
0 . 5000 
0.5200 
0.5400 
0 . 5600 
0.5800 
0 . 6000 
0.6200 
0.6400 
0.6600 
0.6300 
0,7000 
0.7200 
0.7400 
0.7600 
0.7800 
0.8000 
0.8200 
0.8400 
0.8600 
0 • 3300 
0.9000 
0.9200 
0.9400 
B 0 T 1 * 0.9600 

BOTla 0.9800 
6071* 1.0000 

BOT 1 * 1.0200 

BOT la 1.0400 
BOT 1 c 1.0600 
BOT 1 “ 1.0800 

BOT 1 * 1.1000 



TOP assumed l argep 



TOPl> 
T0Pl = 
TOP 1 = 



0 0 7 0 
77pl 
0 0 .7. 3 



than bottom 
BOTl* 

BO T 1 = 

BOT 1 = 

BOT 1 = 

BOT 1 = 

BOT 1 c 
BOT 1 = 

BOT 1 c 
BOTla 
BOT 1 c 
BOT 1 = 

BOTla 
BOT 1 = 

BOTla 
BOTla 
BOTla 
BOTla 
BOTla 
BOTla 
BOTla 
BOT 1 = 

BOTla 
BOTla 
BOTla 
BOTla 
BOTla 
BOTla 
BOTla 
BOTla 
BOTla 
BOTla 

BOTla 
BOTla 
BOTla 
BOTla 
BO” 1 = 

BOTla 
BOTla 
BOTla 
BOTla 
BOTla 
BOTla 

BOTla 

BOTla 
BOTla 
BOTla 
BOTla 
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1000' 
al/5.75 
B/3.3 
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18000 psi 
60" 

STEEL 

4.79" 

1.15" 

18017.24 

39055.43 

39843.44 
22.55' 
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Since sometimes it is desirable to use aluminum alloys 
in the structure of ships, due to their light weight, high 
corrosion resistance and non-magnetic characteristics, in the 
following runs, the material considered is aluminum. The 
necessary change is made in the input to account for the pro- 
perties of this material. 

When considering aluminum as the structural material, 
we must keep in mind that the aluminum hull will be more flexi- 
ble than a corresponding steel ship. The modulus of elasticity 
of aluminum is about one third that of steel, and if scantlings 
were identical, the deflection of an aluminum ship, under a gi- 
ven load, would be three times those of a steel ship. 

However, to obtain a factor of safety on ultimate strength 
equal to that in steel, the area of an aluminum member in direct 
stress, or the section modulus of an aluminum member in bending 
would be about 1.5 times that of the corresponding steel member 
(based on an ultimate strength of about 60000 psi in medium steel 
and of about 40000 psi across the welds of a welded aluminum 
structure). 

A section modulus of 1.5 times that in steel, with the 
same depth being considered in both cases, would result in a 
deflection in aluminum about twice that in steel. However, the 
lower modulus of elasticity of aluminum reduces its ability to 
resist buckling , and for the same factor of safety both the 
section modulus and the moment of inertia in aluminum would in 
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general be somewhat more than 1.5 times those in steel, so 
that the deflection of the aluminum structure would in general 
be somewhat less than twice that of the steel structure. 

The greater deflection would reduce the natural frequency 
of hull vibration, and sowould affect the interaction of slamming 
stress and bending moment stress. It could also affect shaft 
bearing pressures. It has been proposed as an arbitrary basis 
of design that aluminum ships be designed for a deflection about 
1.5 times that of steel ships. This would require a moment of 
inertia of twice the value of that in the steel hull girder, 
and would automatically result in a section modulus substantially 
more than 1.5 times that in steel, and a factor of safety cor- 
respondingly greater than in steel. 
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SUMMARY 



Length 


= 1000' 


Beam 


= AL/ 5.75 


Draft 


= B/3.3 


Depth 


= AL/l4 


DST 


= 18000 psi 


S 


= 120" 


Material 


- ALUMINUM 


TOP 


= 4.27" 


BOT 


= 2.98" 


STH 


= 17958 .79 


SBH 


= 22852.76 


SBCRIT 


= 22768.12 


DITOP 


= 31.43' 


DIBOT 


= 39.99' or 40.00' 
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Looking at the following tahLe we can also observe for the 
aluminum the effect of changing S on the neutral axis location 
and the compressive stresses resulting on thebottom plating. 



CASE 


S 

inches 


TOP 

inches 


BOT 

inches 


SBH 

psi 


SBCRIT 

psi 


DIT0P 

ft 


DIB0T 

ft 


1 


10 


5.35 


0.363 


53000 


49000 


18.09 


53.34 


2 


30 


4.86 


1.00 


41000 


41000 


21.72 


49.71 


3 ; 


60 


4.53 


1.76 


32000 


32000 


25.83 


45.60 


4 


120 


4.27 


2.98 


23000 


23000 


31.43 


40. 00 



(*) This difference may result from inaccuracy in reading 
the plot. 

Here again the values of TOP and BOT may allow for a quick 
evaluation of the weight relationship between the four cases. 

One run was also made in the case of the aluminum for 
AL=500 ft and the results are here included. 
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Intersection not accepted; assumed Top>Bot ended up 
with Bot>Top. 
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SUMMARY 



Length 


= 500 * 


Beam 


= AL /8 


Draft 


= B/3 


Depth 


= AL/9 


DST 


= 18000 psi 


S 


b 

(T\ 

il 


Material 


- ALUMINUM 


TOP 


= 0 . 38 " 


BOT 


= 0 . 6 " 


STH 


= 18089.66 


SBH 


= 14578.14 


SBCRIT 


= 14767.8 5 


DITOP 


= 30.76" 


DIBOT 


= 24.79" 



HOGGING BENDING MOMENT=SAGGING BENDING MOMENT 



DESIGN STRESS IN TENSIOM=DESIGN STRESS IN COMPRESSION 

From the preceding section we can understand the need to 
put some limit on the value of the compressive stresses, and per- 
haps also the inefficiency of considering the value of S (the 
spacing between longitudinal stiffening members) as an input to 
the problem. 

In this section we will consider a bending moment of the 
same magnitude in both hogging and in sagging and the stress 
in either tension or compression will have as upper limit the 
design stress. 

Since the bending moment is the same, it is intuitive that 
the thickness required to resist tension on top and compression 
on bottom (buckling) when the hogging moment is acting will 
also be respectively the thicknesses to resist tension on bottom 
and buckling on top when we consider the sagging moment. We 
would end up then with a thickness on top required by the 
hogging moment TOPH, a thickness on bottom required by the hog- 
ging moment, as well as a thickness on top required by the sagging 
moment TOPS and a thickness on bottom required by the sagging 
moment BOTS. 

Since HM = SM we must obtain 

TOPS = BOTH = tl and BOTS = TOPH = t2 

and we will end up considering the largest tl or t2, to 

be the thickness of the top and bottom simultaneously. 

Assume for instance that tl > t2 then if we decide on tlfor 
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thickness of the top as well as at the bottom, we will 
verify that we have overdesigned because TOPS = tl 
corresponds to BOTS = t2 which plugged into the expression 
for I gives a smaller value than using TOPS = tl with 
BOTH = tl that will give a larger value for I; this larger 
I value when used in 

BM*Y 

STR = 

I 

will give a smaller STR value than the design stress and 
we could use a smaller thickness, may be between tl and 
t2 for both the top and bottom and satisfy the DST requi- 
rement . 

This means that we must incorporate in this analysis the fact 
already known that when BM (hogging) is equal to BM (sagging) 
we will have to consider the neutral axis half way between top 
and bottom which means that the thickness of the top shall be 
equal to the thickness of the bottom for a geometrically 
symmetrical cross section. 

Let's determine the formulas that will best serve this 
case then: 

HM = SM = BM 

DITOP = DIBOT = D*(B*t+D*t)/(2*D*t+B*(2*t)) 
t = TOP = BOT = TS 



I 



2*D*D*t* 



3*b*b+d*d+4*b*d 

D+B 



BM*DITOP BM*DIBOT 



Now we can write 



STR = 



I 



I 



DST 
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and the instability formula 




1-P0I*P0I 



*DST 



E 



Where S is the spacing between longitudinals. 

We see thus that we cannot have more than two equations 
and so we must seek for two variables. A quick inspection and 
one can see that knowing BM, B, D, the material properties and 
the design stress DST we have left to determine what should be 
the thickness of the plating (t) and the corresponding frame 
spacing (S). Solving the equations with this in mind, we obtain 
BM = 1 . 37143*B*T*AL*AL lb*ft 



BM*(D+B) 

DST 

4*D*t ( 3*B*B+D*i>f 4*B*D ) 



which gives 



BM*(D+B) 

t 

4*D*DST* ( 3*B*B+D*D+4*B*D ) 



and also 



BM*Y 



and from 



STR = 



we obtain 



I 



t 



S 



0. 30396* (1-P0I*P0I)*DST 
E 
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These formulas will give for a prescribed set of dimensions 
(ft) and a particular material subject to bending BM (lb*ft) the 
thickness t (inches) and the spacing of longitudinal stiffening 
members S (inches) required. 



HOGGING MOMENT AND SAGGING MOMENT OF 



DIFFERENT MAGNITUDES 

THE DESIGN STRESS IS THE SAME IN TENSION AND 
IN COMPRESSION 



Let's consider now the case where the magnitude of HM 
different than the value for SM. 

We have in this case: 

D* ( B*BOT+D*TS ) 

DITOP = 

2*D*TS+B*(TOP+BOT) 

D*(B*TOP=D*TS) 

DIBOT = 

SAME DENOMINATOR 

3*B*B*TOP*BOT+ ( D*TS ) 2 +2*B*D*TS* ( TOP+BOT ) 

I = 4*D*D* 

SAME DENOMINATOR 



HM*DITOP 

STH 

I 



SM*DITOP 

STS = 

I 



HM*DIBOT 

SBH = 

I 



SM*DIBOT 
SBS = 

I 



BOT = 




TOP = 




S 2 *( l-POI*POI) 



where 



G = 0.30396* 



E 
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Looking at these equations, we know that the result obtained 
for the thicknesses is going to determine the value of DITOP, 

DIBOT and I. If furthermore we restrain the stresses from being 
larger than DST we will undoubtedly want them to end up being 
DST otherwise if they are smaller we are not using the material 
efficiently. 

We thus have 

DST = STH = HM*K1 DST = STS = SM*K1 where K1 = DITOP/l 

DST = SBH = HM*K2 DST = SBS = SM*K2 where K2 = DIBOT/l 

This tells us that if for instance HM > SM 
we should be considering DST = STH = HM*K1 

and DST = SBH = HM*K2 

because the other two equations will determine lower stresses 
and are thus within the requirement of not exceeding the 
design stress, being automatically satisfied. 

Now from these equations we can see that if the design stress 
required for the top, (in tension for the case of HM>SM; note 
that for the case SM > HM the other two equations should be con- 
sidered and the conclusions are similar), is going to be the 
same one as required for the bottom ( in compression ) , we must 
have K1 = K2 and so DITOP = DIBOT and the neutral axis will 
again end up being located halfway between top and bottom for 
a geometrically symmetrical cross section. 

We can here use the same two equations derived in the 
previous section to determine (t) and (S). The moment to be 
considered here will be the larger of HM or SM instead of BM. 
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HOGGING MOMENT LARGER THAN SAGGING MOMENT 
DESIGN STRESS IN COMPRESSION LARGER THAN DESIGN STRESS IN TENSION 

Thesituation we are analysing in this section 

HM > SM 
DSC > DST 

is perhaps more interesting than the previous ones because its 
study provided us with a much simpler method for determining the 
thicknesses and the neutral axis location; we call here the atten- 
tion of the reader to fully understand the steps of this section 
because the following sections till the end of this work will use 
this same method of analysis applied totheir own particular case. 

Indicated in the next diagram are the ranges within which we 
may stress the material in order to satisfy the requirement of 



exceeding DSC or 


DST. 

DSC 


r* DST— | 




1 


i 

) 

t 

















COMPRESSION TENSION 





Let's consider now that we will try to determine the best 
position for the neutral axis by first gaining some insight into 
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the geometry of the stress distributions and the stress limits 



already defined. 



Unacceptable 




Prom the previous diagram we may easily understand that 
if we decide to stress the material to its maximum in tension 
and compression when under the effect of a sagging bending 
moment (SBS - STS) we end up with a neutral axis closer to the 
bottom due to the fact that DST < DSC and we see that when consi- 
dering the hogging moment the only possible alternative left 
is to limit the stress caused by this moment in tension on top 
STH=DST; in this case the compressive stress on the bottom SBH 
would automatically fall inside the acceptable range. 

If for the hogging moment we choose to allow the stress 
in compression on the bottom to attain its maximum value SBH=DSC 
we see from thediagram that we are forcing the STH to fall out- 
side of the acceptable range. We are left in this situation with 
an acceptable solution so far and we. must investigate it further. 
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At this stage we have 




We may see from the diagram that we are asking for SBS=STH= 
=DST to have the same value and at the same time to have the 
neutral axis located closer to the bottom DITOP> DIBOT. 

We note that the equations that give STH and SBS are 

STH = HM*DITOP/l 
and SBS = SM*DIB0T/I 

and since I is the same STH can never be of the same magnitude 
as SBS because 

STH = HM*DIT0P/I > SBS = SM*DIBOT/l 

both because DITOP is being set larger than DIBOT and also 
because we are considering HM > SM to begin with. Thus we are 
left with the alternative of making DIBOT larger than DITOP 
to make it possible to satisfy 

STH = HM*DIT0P/I = SBS= SM*DIBOT/l 
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Thesituation would then be 




Besides this is only natural because what we are doing here 
is to make the material to be stressed to its allowable limit 
while resisting the larger of the two moments; hogging in this 
case. This situation will be worked out here as an example and 
it seems to be quite clear that similarly, cases where SM> HM or 
DST > DSC or other situation like DST^DSC with HM > SM or 
DST<DSC with SM > HM could be solved applying a similar reasoning. 
So, having HM>SM and DSC > DST we should consider the smaller 
design stress ? in this case DST ? and have the material stressed 
to this upper limit. Why?, because if the larger is considered, 

DSC in this case, there is no way we can have the material stre- 
ssed to this higher limit in both hogging and sagging while con- 
sidering the smaller we may attempt to do so. This can be better 
understood from the two next diagrams. 
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So, considering the lower design stress, in this case DST, 
we will stress the material to this maximum in both hogging and 
sagging thus giving 

STH = HM*DIT0P/I = DST 
SBS = SM*DIB0T/I = DST 

Note here that we already know that the neutral axis will 



be closer to the top DITOP < DIBOT which means that we must obtain 
a thickness for the plating at the top larger than at the bottom and 
so TS=TOP which means that we will be using the following for- 
mulas 



DITOP = 



D*(B*BOT+D*TOP) 

P 



D*TOP*(B+D) 

DIBOT = 

P 



I 



4*D 2 *T0P* ( ( 3*B 2 +2*B*D ) *BOT+ ( D 2 +2*B*D ) *TOP ) 



P 



where P =(2*D+B)*T0P+B*B0T 
we have also 

STH = DST = HM*DITOP/l 
SBS = DST = SM*DIB0T/I 
SBH = DSC = HM*DIB0T/I 

From these three equations we may combine the two first to obtain 
HM*DITOP = SM*DIBOT 

and combining the second and third we obtain 
DST*HM = DSC*SM 

After combining these two equations we obtain 
DST*DIBOT = DSC*DITOP 

This equation could have been obtained more easily just 
by looking at the diagram we are considering and notice the 
similar triangles 
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DST 




Note also that what we obtain from this equation is to 
determine the N.A. location by making DITOP and DIBOT proportio- 
nal to DST And DSC. 

Proceeding with this equation we obtain: 

DST*(D*TOP*(B+D))/P + DSC*(D*B*BOT+D 2 *TOP)/P 
TOP*(DST*B+DST*D-DSC*D) = DSC*B*BOT 
and making 

AC = DST/DSC+D*(DST-DSC)/(DSC*B) 
we obtain 
BOT = AC*TOP 



We need now an equation to stress the material to its design 
limit and this can be done by either using 
STH = DST = HM*DIT0P/I 
or SBH = DSC = HM*DIBOT/l 

Prom these two equations I prefer to use the second especially 
because it will simplify the calculations that follow as can 
be understood from the simpler expression that gives DIBOT as 
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compared with the one that gives DITOP 
We will then consider here 

SBH = DSC = HM*DIB0T/I 
which gives 

HM*D*TOP*(B+D) 



DSC = 


4*D 2 *T0P ( ( 3*B 2 +2*B*D ) *BOT+ ( D 2 +2*B*D ) *TOP ) 



Simplifying we obtain 



TOP = 


HM* ( B+D) 3*B 2 +2*B*D 

= *bot 

DSC*4*D* ( D 2 +2*B*D ) D^+2*B*D 


or making 





A = HM*(B+D)/(DSC*4-“11^D 2 +2*B*D) ) 



AB = 


( 3*B 2 +2*B*D ) / ( D 2 +2*B*D ) 



we have 



TOP = 


= A-AB*BOT 



replacing here 



BOT = 


= AC*TOP 



we obtain 



TOP = 


A 

1+AB^AC 



Having already found TOP and BOT we may use BryanS formula to 
determine S; since the thickness of plating is for this case 



smaller at 


the bottom we must consider this thickness and 


along with 


it DSC because this will be the larger compressive 


stress the 


bottom will be subjected to. 
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From 



DSC 



E 

o. 30396* (l-POI*POI) 



2 

BOT 




We obtain 

1 

S = ■■ - ■ *BOT 

l-POI*POI 

\ 0. 30396* *DSC 

N E 

or making 

1 

Z = 

1 -P 0 I*P 0 I 

\ 0.30396* *DSC 

N E 

we have 



S = Z*B 0 T 

A computer program was written and applied to two parti- 
cular sets of dimensions and two different materials namely steel 
and aluminum. The results are presented next in the computer 
printouts . 
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HOGGING MOMENT AND SAGGING MOMENT OF EQUAL MAGNITUDE 
DESIGN STRESS IN TENSION AND DESIGN STRESS IN COMPRESSION OF 

EQUAL VALUE; LATERAL LOAD ALSO CONSIDERED 

The lateral load being here considered will be a "head" 
of salt water of eight feet on the top plating (HT0P=8 ft) 
and a "head" of salt water equal to the draft on the bottom 
plating ( HBOT=DRA ) . 

We will use here the simplified approach as suggested 
in the following note; this simplified approach is at least 
conservative . 

Note: Bleich, from his experiments and theoretical studies 
on the effect of normal pressure in the buckling strength of 
thin panels, concluded that the buckling strength is increased 
when combined with a sufficiently large normal load, and that 
the increase is greater for thes imply supported edges situation 
than with fixed edges. For heavier plating , where thedeflection 
does not exceed half of the plating thickness, the increase in 
buckling strength resulting from normal load is negligible. 

It was therefore suggested, for simplicity, that the critical 
stress be computed for in plane loading only. 

The tertiary stresses we are interested in considering 
here are those exerted at points B directed longitudinally ; these 
will combine with the longitudinal bending stresses or primary 
stresses as described in a previous section. 

Since now we have DST=DSC and HM=SM=BM we may think of 
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a diagram for the stress distributions as indicated next figure: 




Prom this diagram we may understand that the neutral axis 
should be located closer to the bottom which means that the thi 
ckness at the bottom will be larger than at the top BOT>TOP . 
So, we may already write the following formulas where TS was 
replaced by BOT 

DITOP = D*(B+D)*BOT/P 

DIBOT = D*(B*TOP+D*BOT)/P 

I = 4*D*D*( ( 3*b^+2*b*d)*top+(d2+2*b*d)*bot)*bot/p 



where 



P = 2*D*B0T+B*(T0P+B0T) 
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We may now think of the following relationships to deter- 



mine TOP and BOT. 



STH 




SBH 

STH*DIBOT = SBH*DITOP 
or knowing that 

S ? 

STR3T = 0.152222*HTOP*( ) d 

TOP 

S ? 

STR3B = 0.152222*HBOT*( ) d 

BOT 

and anticipating 
S = Z*TOP 

where 



1-P0I*P0I 




*DST ) 



we have 



2 



S 



Z 



TOP 



S 




If we now make ( 



BOT 



in order to simplify the equations, we see that when we do this, 
since 
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■) 2 < (• 



BOT 



TOP 



we are increasing the value of STR3B and. this means that since 
DST = SBH+STR3B 

we are decreasing"slightly" the value of SBH that could be allo- 
wed to reach a larger value; so , the substitution 




BOT 



is on the safe side and has the advantage of simplifying the 
equations considerably. 

So we may write 
STH = DST-STR3T 
SBH = DST-STR3B 
STR3T = 0.152222*HTOP*Z 2 
STR3B = 0.1 5222 2*HB0T*Z 2 
A = DST-STR3T 
AB = DST-STR3B 
so 

A*DIB0T = AB*DIT0P 

and we obtain replacing DITOP and DIBOT by their equations, 
and simplifying 

AB*(B+D)-A*D 

TOP *B0T 

A*B 



or making 
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AC = (AB*(B+D)-A*D)/(A*B) 
we finally obtain 
TOP = AC*BOT 

Now in order to stress the material to its design stress 
we may either consider STH, STS, SBH, SBS. 

Due to ease of calculation we may either choose STS or 
STH; the expressions will be exactly the same, I will consider 
here 

STH = BM*DITOP/l 

and 

STH = DST-STR3T = A 

Substituting for DITOP and I we obtain 



BM*D*(B+D)*BOT 



4*D 2 *B0T* ( ( 3*B 2 +2*B*D ) *TOP+ ( D 2 +2*B*D ) *BOT ) 



and after simplifying we arrive at 



BM*(B+D) 



3*B 2 +2*B*D 



BOT 





*TOP 



Making 



AD = BM*(B+D)/(A*4*D*(D 2 +2*B*D) ) 
AE = (3*B 2 +2*B*D)/(D 2 +2*B*D) 



we obtain 



BOT = AD-AE*TOP 
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substituting 

TOP = AC*BOT 
we obtain 



AD 

bot 

1+AC*AE 



Considering now the instability formula (BRYAN) we notice 
that since the top will have the thinner plating we will consi- 
der TOP stressed to the maximum allowed DST. 

So we start with 



Z = l/( 



N 



1-P0I*P0I 

0.30396* *DST ) 



S = Z*TOP 



The total stresses will then be given by 
TSTH = BM*DIT0P/I+STR3T 
TSTS = BM*DIT0P/I+STR3T 
TSBH = BM*DIBOT/I+0. 152222*HB0T*S 2 /B0T 2 
TSBS = BM*DIB0T/l+0. 152222*HB0T*S 2 /B0T 2 



The computer printouts for the same cases analysed in theprevious 
section are presented next. 
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C *al* IS the LFNGHT 

AL»5P0 • 

R'AL/8 • 

C *DR A * IS THE oRAFT 

DRAt=B/3. 

0~ AL/9 • 

C *POI* IS THE oolSSON's RAT I° OF THE MATERIAL 

C FOR STEEL 

pUl«0. 3 

C *E * IS THE YOUNG'S MOpULUS OF T H F- MATERIAL 

C FOP steel 

E*30 ♦ * 1 P • **6 

C *DST* IS THE nESlGN SyRESS 

Oi=T = E5000 . 

c *pm* is the bending moment 

B M =1 .37143*B*prA*A l.*A|_ 

W*0»303B6M1‘-dOI*PoI>*DsT/E 
U = SORT ( l' 1 ) 

Z=1 ./u 

♦htop* IS THE HEAD oE wajep on deck 
H TOP s 8 • 

*HSOT* IS the head Of Water ON BOTTOM 
hbot=dra 

STP3T-0* 15222?* WTOP*Z.Z 
STR3B = 0 • 1 52?2?* u R0t*Z*Z 
A«DST-STR3T 
AB«0ST-STR3B 
AC»(AB+(R+D)- a *D)/(A*r) 

AD«sPM* ( B + D ) / ( a*4«+D*( 0*D + 2«*B*D ) 1 
aEs(3.*B*B+2.*r*0»/(D*D+2.*B*0) 

C *TOP* IS THE tHICKNFSq OF THE tOP PLATING 

C *90T* IS THE THICKNESS OF THE BOTTOM PLATING 

BOTbAD/ < 1 . +AC*AE ) 

TOPbAC*BOT 
WRITE ( 5/ 10 ) TOp 

If FORMAT; « 1 * 'TOPs'plO. 4) 

WRITE(5>20)BOt 

20 FORMAT!' ' / 27 yi ' BOTc ' cl0 . 4 ) 

C *S* IS THE SPACING OF LONGITUDINAL STIFFENING MEMBERS 

S*Z*TOP 
WR I TE ( 5> 30 ) S 

30 FORMAT;' ' ,2?y, »S= 'Fig. 4 ) 
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W r ITE(5>40)DIto p 
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C bDIBOT* IS THp DISt a N c E OF THE BOTTOM FROM THE NEUTRAL AxIS 

DIROToD*(B*Top+0«ROT)/P 
WR!TE(5i50)OlqoT 

50 FORMAT;' ' , 1 8y , ' D I ROT, ' F 1 0 . 4 ) 

C ^M I * IS THE MOMENT Of INEPTIA OF THE X SECT I OK ABOUT N#A. 
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TOP= 0.1671 

BOT= 0.8?83 

S= 11.0049 
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Note here that because we have in this run BOT>> TOP the 

approximation made to simplify the formulas, namely 
S 

Z =- 

BOT 

becomes very inaccurate compared to what the actual value 
should be 

S 

Z 

TOP 

this is the cause for obtaining such low stresses at the bottom. 
We might refine the result for this case as is done for the 
"1000 ft steel case" also included in this section which presen- 
ted still another difficulty. Since all steps to have the bottom 
stressed to DST are explained in detail for such case they will 
be omited here. 
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This case gave a result close to the correct solution with 
the bottom stressed to almost DST. 

Next is presented the "1000 ft steel case" where still 
another difficulty was encountered and a procedure to have the 
bottom stressed to DST is exemplified. 



171 



STRENGhT Of SHIPS-ON the NEUTRAL axis location page 



•04 


C 




* A l * IS THE l-ENGHT 


■04 






AL=10O0« 


'K 






p-AL/5 . 75 


as 


C 




* DR A * IS T^F DRAFT 


18 






DR A = D/3 • 3 


Ik 






D-AL/14* 


30 


c 




* P 0 1 * 10 T.-.E POISSON'S RATIO OF THE MATERIAL 


30 


c 




for stell 


30 






pO I =0 • 3 


38 


c 




*B* IS The YOUNG'S MODULUS OF THE MATERIAL 


38 


c 




FOR STEEL 


38 






E *"30 * * 1 0 • * ♦ 6 


‘•C 


c 




* C S T * IS Tme DESIGN stress 


kC 






D ? T =25070 • 


5k 


c 




♦Eh* IS THE BENDING MOMFNT 


bk 






BE =1 •37l43»B*DRA*AL*AL 


6C 






W=C.3K396*(l.-P0I*PCI)*DST/t 


eC 






U = S 0 R T ( V» ) 


98 






z = l */U 


Ak 






* H T OP » IS THE HEAD OF WATER ON DECK 


Ak 






HT 0P = 8 • 


4 C 






*HBOT* IS THE HEAD OF WATER ON BOTTOM 


AC 






uBoT =DRA 


ck 






S T R 3 T = 0 • l5?2?2*HT0P«Z*Z 


C8 






STR3Bs0. l522?2*h?0T*Z»Z 


JC 






A=DST-STR3T 


1 8 






AE=DST -STK3B 


'k 






AC=(AB*«f5 + C)-A»D)/(A*R) 


1C 






AD=BM*(B+D)/(A*4.*D+(D*D+2.*B*D)) 








AE=(3»*B*6-t-2.*F*p)/(D*D + 2»*B*D) 


Ik 


c 




* T 0? * IS The THICKNESS OF THE TOP PLATING 


Ik 


c 




* B 0 T * IS T-E THICKNESS OF THE BOTTOM PlaTING 


Ik 






BOT = AD/ ( 1 • + AC* AE ) 


)c 






T0P=AC*60T 


. 

:8 






write < 5> i0 j top 


!k 




10 


FORMAT ( ' 1 ' 0 2tfX, ' TOP= ' F 10 . 4 ) 


IE 






WRITE(5*20)BOT 


lA 




20 


FORMAT ( ' ' ,2/X, 'BOT='F10.4 ) 


i>4 


c 




* S * IS THE SPACING C-F LONGITUDINAL STIFFENING MEMBERS 


i)k 






S=Z*TOP 


l>0 






WPITE(5>30)S 


rc 




30 


FORMAT!' S 2?X, ' S= ' F 10. 4 ) 


Dk 






pe 2 »*D*BcT + B*( TOP+OOT ) 


*8 


c 




- DITOP* IS ThF. DISTANCE OF THE TOP FROM THE NEUTRAL AXIS 


J >8 






DITOP=D*(B«BoT+D*RCT)/P 


l)C 






P I T E (5/40 ) D I T 0 P 


1 '8 




40 


FORMAT!' ' > 18X, 'DTTOP='F10.4 ) 


Ik 


c 




*D I BOT * IS ThE DISTANCE OF THE BOTTOM FROM THE NEUTRAL AXI 


1 k 






DlBOT=O*(B»T0P+D*GOT )/P 


1 18 






WRITE ( 5/50 ) DIBOT 


2 *4 




50 


FORMAT ( ' ' , 18*> ' DIBOT= ' F 10 . 4 ) 


*'0 


c 




*GM I * IS The moment of inertia OP THE X SECTION ABOUT N • A • 



STRENiiHT Op-' SHIPS-ON ThE NEUTRAL AXIS LOCATION 



172 

PAGE 2 



|70 

i£0 

.E0 

U 

E? 

£0 

£0 

£0 

FA 

: 8 

6A 

7E 

9a 

5a 

DP 

EC 

28 

6 A 

AA 

C0 

Ok 

'6 

10 



C 

c 

c 

c 

c 

c 



60 



70 



80 



90 



00 



GMI=A.*D*D*( |3.*P*B+2«*B*D)*TOP 
-TSTh*IS The STRESS ON TOP OUR 
*TSTS»IS T - E STRESS On TOP DUE 
» T S E3 H * IS Ti-E STRESS ON POT DUE 
».TSBS*jS The STRESS ON POT DUE 
■^STCPIT* IS THE CRITICAL STRESS 
-SRCRIT+ IS THE CRITICAL STRESS 
TSTH = E3M*Dl tOp/GMI-*-STR?T 
TSTS = tJM*DlTCR/G^I+STR3T 
STCRIT=TOP»ToP*E/(0.3C3y6*S+S* ( 
WE T TE ( 5/60 ) TRTu 

FORMAT ( ' ' , 19*/ ' TSTH= I F 10*3 ) 

WRITE ( 5/ 70' ) TbTS 

FORMAT)' 1 / 1 g X / , TST?='El0«3) 

WRITE(5/80)STCRIT 

FORMAT) ’ '/17X/»STCRIT='F10.3) 

TS8H = PM-»DI GOT/GMI+C . 1 52E22* HfiOT 

TSUS = SM*DI HOT/GM I+.-> . 1522 22 ‘HOOT 

SB CR I T =BQl *B0T*E/( 0.30396 -S*S* ( 

wRITE ( 5/ SO ) T SBh 

FORMAT)' ' . 1°X/ ' TSEu= ( F10.3 ) 

wR I T E ( 5 / 1 0 , ' ) T S R S 



*ROT + (D*r:+2.*P»-P)*B0T*B0T )/P 
TO HOGGING 
TO SAGGING 
TO HOGGING 
TO SAGGING 
ON TOP 
ON B G T 



1 . -PG I *P0 I ) ) 



*S*S/ ( 6 0 T * B 0 T ) 
*S*S/ ( 60T*B0T ) 
1 . -P0I*P0I ) ) 



FORMAT ( • ' , 1?X , • TSBS- 

wRI TE ( 5/ 1 10) SBCRI T 



' P10.3 ) 



t?C 


: 10 


format ( 


' ' / 1 7 X / 


•SBCRIT= 


' F 1 C • 3 ) 














END 
















(S3 


• u 


26AC CV] 


AL 


065 A CV3 


B 


0 6 5 C CV] 


DRa 


06 6 A C V ] 


D 


CV3 


POI 


06 7 A [ V 3 


E 


0000 CS3 


• R 


0 6 8 A [V] 


DST 


;)68C CV] 


gn 


UVj 


r* 


06 A 0 [ V 3 


U 


0/00 CS3 


SORT 


06 A A CV] 


1 


06 A 8 CV] 


rlTOP 


CV3 


HBOT 


06BA [Vj 


STR3T 


06«C C V J 


STR3B 


0 6 C U CV] 


A 


O' 

0 

4m 

r*-> 

< 

1 • 


Ab 


(Vj 


AC 


06DC c v : 


AD 


06Pf» CV 3 


AE 


0 6 F 8 CV3 


BOT 


06P c : v ] 


TOP 




10 


0000 [S3 


(il 


023 A CL 3 


2 0 


0 7 v 0 C V ] 


S 


027 C CL3 


3 A 


HVj 


P 


0708 C V 3 


Dl TOP 


02R0 CL 3 


A 0 


0 7? C [V] 


01 BOT 


0.35 A CL 3 


5/ 


Cv3 


GMI 


■0 7 1 C C V 3 


tsth 


07?-7 CV3 


TSTS 


0 7 2 A CV] 


STCPIT 


U A 6 A C _ 3 


6 v) 


s:l3 


70 


0 A p0 [L 3 


80 


072? Cv3 


TSBH 


072C CV] 


TSBS 


0731.- C V] 


SB C* I 


^•L j 


90 


05F6 C L 3 


100 


062C Ci_3 


110 


0000 lS] 


. V 







:o 



’ >R Af- LAPELS: 



70 


* M A I U * 


2E 3 A 


. V 


2C6A 


.COMP 


298 C 


ALOG 


2EA8 


a i 


2 7 A A 


. R 


'32 


• ZERO 


2C90 


$6 


27EA 


. A 


2D8E 


•M£S 


2C0C 


. w 


2kS h 


EX< 


9 A 


AINT 


2896 


SORT 


2CD2 


.RARG 


2D2E 


.5 


2D0A 


s8 


2036 


• E 


36 


• 0 


2E0C 


.U 


3F 6 A 
















Y- 


points : 






















AA 


• R 


2 7E A 


. A 


2896 


SORT 


298C 


ALOG 


2A9A 


EXP 


2B5 A 


IP 


3C 


• w 


2C6A 


. Comp 


2C9P 


$6 


2CD2 


. R A p G 


2 00 A 


«8 


21. 2 E 


.5 


32 


• ZE^O 


2036 


•ERCnT 


2038 


• 0 


2D«E 


. me s 


2E18 


• U 


2E3E 


. V 
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T OP = 
BOT = 
S = 

OlTOP= 
0 1 BOT = 
TSTH = 
TST3 = 
S T CR I T = 
TSBH = 
TSB3 C 
SBCRI r= 
ND 



VIC OPERATING SYS 



6 • 1 7 4 8 
- 5 • 5 5 ? 1 
40 * . 7124 
1 4 2 . 0 6 1 4 

-70.6328 
249S9 . 99? 
24999. S92 
24999.983 
3324 1.342 
33241 . 84 / 
2 2 1 3 . 25 / 



version 1 REVISION 



012 P3/04/74 GENERATED 



)B 


HANDLING CHARGE 


* . 3b 


/ jnB 


.35 


25 


LINES PRINTED FR2 


* 1.25 


/ K ln 


•16 


5 4 


Cards read 


s 1.50 


/ * CD 


.13 


*0 


plotter vectors 


s .25 


/ 10 0 7 


. 0 0 


L- 1 


model 70 SECONDS 


*25.00 


/ HOUR 


* 1 5 


'■2 


v ODEl 30 seconds 


*12.50 


/ n r u r 


. 0/ 






total 


CHARGE 


? .79 


N IK 


452 14731 LOGGED 


CUT 0 


4/30/74 


21126. * 27.07 LEFT AFTER 



04/30/74 



53 LOG I J3 * 
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It does not come as a surprise that the program does not 
work for the "1000 ft steel case" since in this situation we 
obtain STR3B DST which cannot be accepted and also this 
program is not prepared to handle. 

AL = 1000 

B = AL/5.75 = 173.913 
DRA = B/3.3 = 52.7 
poi = 0.3 

E = 30*10 6 
DST = 25000 

BM = 1.37143*B*DRA#AL 2 = 12569. 45*1 0 6 
w = o. 30396 *(i-poi 2 )*dst/e = 230 . 503 * 1 o ~ 6 

U = W = 15.182*10 “ 3 

Z = l/U = 65.867 
HB0T = DRA = 52.7 

STR3B = 0. 152222*HB0T*Z 2 = 34803.57 DST 
Since we understand that the program should work for DST^STR3B 
we make the next run with DST = 35000 to prove that in such case 
a solution is obtained^and after proceed to obtain the actual 
solution for the case we were considering with DST = 25000. 
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page 



1 



•04 


C 




* A L * IS TNfr L Ef,GHT 


•04 






A L = 1 H «' •<* • 


•0C 






B = l L / b • 7 5 


18 


c 




* DP A * IS T...K DRAFT 


18 






0 K/ =B/3 • 3 


24 






D = : L / 1 4 • 


30 


c 




«PPI* IS T-E POISSO' '5 PATIO OF THE MATERIAL 


30 


c 




P OF STEEL 


3(7! 






P 0 1 = P • 3 


3s 


c 




* F » IS THE YOUNG'S • ODULUS OF THE MATERIAL 


38 


r 

•w 




POP STEEL 


38 






* 10. **6 


K 


r 




* D S' T * IS T^E Dp 3 I ON STRESS 


*C 






0 s 1 = 3 5 0 0 • 


54 


c 




*B‘i* IS T mf BE* DING MOMENT 


54 






8F = 1 • 3 7 1 4 3 * B » D R A # AL *AL 


bC 






Y =C. • 30 396 * 1 1 . - P’O! I ) *DST/E 


SC 






U»ffJ»T ( W ) 


98 






Z=) • /D 


A4 






* H 1 OP* IS IMF HEAD "F WATER ON DECK 


■A 4 






H T(’P = 8 • 


K 






»MPOr* IS TP'E heap of WATER on bottom 


AC 






H E P T = t ; R A 


o4 






STP?T = t . l!>2?22*MTpP*Z*Z 


C* 






ST D 3B«0 • 15SB22-MB0T- Z»Z 


3C 






A = r >ST - C T R 3 T 


'•'8 






A B = D S T -S TP pP 


” 4 






AC= ( AT* ( B+ '■ ) - A * D ) / ( a*D ) 


■2C 






A D = B M * ( F + 1 ' ) / ( A * 4 * * 0 ^ ( D * D + 2 • * U * D ) ) 
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aE=(3«*B*F+c’* !,c F*1 )/ 0-'D-*-2» *B*D ) 


*4 
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• TOP* IS 1 ‘E THICKNESS OF THE TOP RLATjrG 


1:4 


c 




* BO T * IS T„e THICKNESS OF THE BOTTOM Pl-' TING 


1:4 






P 0 T = A D / ( 1 • * A r. * A E ) 


»>c 






top=ac* r o t 


t: 8 






WRITE (5> 10 ) Top 


e:4 




10 


FORMAT ( • 1 1 JC?X, ' T C P = ' F' 1 0 • 4 ) 


r 






WR tTF. ( 5 > 20 1 BOT 


14 




20 


FORM A T ( 1 ' , 2?X, • rOT= ' F10 • 4 ) 


' ’4 


c 




« S * IS the SPACING of LONGITUDINAL STIFFENING members 


4 






S = Z * T G P 


5.0 






WRITF ( Si 30 )S 


*C 




30 


pORMAT(' ' , 22X, • Ss • P10.4 ) 


i'4 






P-2» M D*BOT + B*(TOP*-BCT ) 


•>8 


c 




*OI T O c * is THE DISTANCE op the top from the neutral AXIS 


« 8 






DITOP = DmB*BoT + D*POT)/ d 


e 'C 






write ( 5/ 4 0 )DlTc p 


t 8 




40 


FORMAT ( ' ' , 1 RX , • 01 TOP= ' F 1 0 . 4 ) 


« 4 


c 




* D I BO T * IS ThE DISTANCE OP THE BOTTOM FROM THE NEUTRAL AXIS 


* 4 






DlBOT=D* ( B-i c p +D»PC T )/P 


* s 






WRITE! 5/53 iDlBCT 


* 4 




50 


FORMAT ( ' ' , ivX , ' DIBOT= ' F 1 0 . 4 ) 


* 0 


c 
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pace 
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GMI=4.*D*Dv((3.*p*P 42.*B*r)*T0P*B0T + (D»O + ?**B*D)*L 1 Pr*B0'I)/P 


E 0 
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♦ tsth* IS 


T t. E STRESS ON TOP DUE TO HOGG I \ G 


E 0 


c 




* T S T S * IS 


Tv-E S t RES-R ON TOP DUE T() SAGGt‘-G 


E 0 


c 




* T SRH* IS 


Tuf s^Rffc; piy pOT DUE To HOGGi-G 


£0 


c 




•TSBS* IS 


Tuf STRESS ON 'BQT DUE T Q SAQGt>G 


E 0 


c 




*STCF IT* 


I*. THE CRITICAL STRESS ON T0 P 


E 0 


c 




*5BCR I T * 
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E 0 






t S T M = f' M * C I ■> 0 P / r j M I + j; ' R 3 T 








TSTS = B;i*rl T 0p/C-MI+S’R3T 
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STCRT1 = T 0 P* T(.'P*E / ( 0 .3(?396«S*S* ( 1 .-POI*POI ) ) 
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WRITE ( b; 60 ) TeTP 
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60 


format ( • 


* , 1°X , * TS T ( j= i F 10 . 3 ) 


7E 






WRITE ( S / 70 ) T STS 


9A 




70 


FORMAT ( • 


' . 1 9 X , • T c t c = ' F 1 0 . 3 ) 
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WRITE <5, 


60 1 STCR I T 


70 




80 


FORMAT ( • 


' , 1 ?X, *STCRIT«'F1P.3) 


kc 






TSB u = f : M* 


DIpOT /GMJ. +0. 1S2???*HB0T*S«-S/ ( B0T*B0T ) 


00 

rvi 






T S B S = c M * 


DI'-OT /GMJ 4-w . i 5222?* U0O T *S* S/ ( brT* GOT ) 








SBCRIT =B 


O r *POT*E/iy.30 396*S*S*( 1 . - PO I * PC I > ) 
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f?C 
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FORMAT ( ' 
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AL 
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06 a 0 c v : 
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AC 
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AD 
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10 
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VI 


P 
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0 2 F c CLI 


VI 


GMI 
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ISth 


77 2 ' C v I 


3 U 


70 


0400 CL I 
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07RF Cv I 


c LI 


90 
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062C C L I 



3 L 



' Fie. 3 ) 



B 


O 6 5 C CV] 


OR A 


0 6 6 4 C V ] 


D 


. R 


0684 CV] 


DST 


V 6 * C C V ] 


PN 


SORT 
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7 . 


0 6 A 5 C V ] 
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STR 3B 


06C0 CV] 


A 


06 c<* cv: 


/ p 


AF 


06FP CV] 


BOT 


0 6 f c c v : 


T C ° 


20 


0700 CX] 


S 


02 7 C CLI 


3? 


40 


070C CV] 


DIBOT 


0 3b '4 CLI 


P7 


TSTS 


0724 CV] 


STCR IT 


04 64 cl: 


f 7 


TSBH 


072C CV] 


TSBS 


0 7 3 " C V : 


S P C R I T 


1 10 


000 vl rS] 
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•MAIN* 


see a 


• V 


2CE A .COMP 


B2 


•ZERO 


2D 10 


*6 


27R4 .A 


1A 


A I NT 


2896 


SORT 


2052 «RAPG 


84 


*8 


2DB6 


• C 


2E8C .U 



►y- 

A 4 


points : 
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. V 
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•5 
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2 A 94 
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2 A 9 4 


E X c 


2D 1 0 


£ 6 


2052 


. R A p G 


208 4 


$ 3 


2DB8 


• C 


2E0E 


• MES 
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. U 
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N'D 



TCP* R.4626 
BCT* 9.7569 
S - ? 5 , 7 5 0 9 

D I T C P = 5? . ni«6 

D I BC T- 57,5090 
TST H = 3*. 9?9 . 977 
T S T 5 = 3 (i 999 • 977 
STCRIT= 3<>C/e./t6 
T S f- H = 1 1 c 6 • 5 4 7 
T S E 3 = 1 196 . 5^7 

Our^nx- 



Vic OPERATING 5 Y f T f * V F: R S I 0 N 1 REVISION g 1 2 03/04/74 O £ N t R A T F l 



3B HANDLING CHARGE 
25 LINE? PRINTED PR? 
S<t C / R :vs READ 

i?. Pi OTTPR VICTORS 
15 V CDEE 70 SE CO N r 9 
?? k 'CDEL F 2 SECONDS 



? .35 


/ 


J-B 


» 35 


* 1.25 


/ 


K LN 


*16 


^ 1.5;’ 


/ 


< CD 


•13 


s .25 


/ 


1 */^0 


• 00 


frf).?'' 


/ 


Mn UR 


• ie 


s 1 2 . 5 ,* 


/ 




.00 


TOTAL 


CHARGE 9 


•74 



IR 490 14731 LOGGED "UT 04/P9/74 20:46. S 32,4s LEFT AFTER 



7> 4 / ? 7 / 7 4 



4 ? LOOPS. 
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Now let's consider our situation with DST = 25000. Once 

the situation is such that due to the approximation made in 
S 

Z 

BOT 

we end up with a STR3B>DST we must consider a value Zl lower 

than Z to bring STR3B down. This can be also understood since 

it is acceptable to make here 
S 

= Z 

TOP 

but because we obtain with the program as is T0P«B0T, making 

the approximation 
S 

= Z 

BOT 

makes Z attain a much larger value than it actually should. 

We must then consider Zl< Z = 65.8 and may start by trying 
Zl = 35 

The computer printout in this case is as follows. 
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3'’ 




Pi I =f .1 
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This time we obtained, a realistic but still unacceptable 
result; we notice that the bottom is being stressed to 31000 ps 



when it can only go up to 25000 psi; besides we assumed 
S 

zi 35 

BOT 

and end up with 

S 225.4 

= 44.46 

BOT 5.07 

This means that the Zl value should be higher. 

We may now try Zl = 40; the computer results for this case are 
shown next. 
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In this case we assumed. Z1 = 40 and. ended up with 

S 199.5 

= = 15.09 

BOT 6.59 

and a stress at the bottom of only 19000 psi which is too small. 



We may now try to use a better Z1 value and can do the 
following: 





^ The results obtained with this Z1 value are the following: 
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The results show that we have obtained. 

S 213.9701 

Z1 ■ = 37.265 

BOT 5.7^18 
while it was assumed 

Z1 = 37.617 

We could further correct the result until Z1 assumed and Z1 actual 
would coincide but it seems that there is no need to refine more 
the solution for the purpose here involved. 

Note that this procedure should also be followed for the other 
cases in order to bring the stress at the bottom to its maximum 
allowed DST for such cases as well. 
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This 1000 ft aluminum case shows again a stress 
bottom much lower than DST. It should be corrected as 
"1000 ft steel" in the previous situation. 



at the 
done for 
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HOGGING MOMENT LARGER THAN SAGGING MOMENT 
DESIGN STRESS IN COMPRESSION LARGER THAN DESIGN STRESS IN TENSION 
LATERAL LOAD ALSO CONSIDERED 

With the help of what was understood for the similar si- 
tuation without the lateral load being included, we may here 
proceed at once to look at what the diagram of stress distribu- 
tions might look like. Actually, the diagram that will be con- 
sidered first is basically the same with only a small change 
due to the addition of the tertiary stresses; since these will 
probably be of the same magnitude on both top or bottom and 
assuming that their actual value does not differ considerably, 
what we may start looking into is a diagram as follows: 




We would consider here the same basic distribution of longitu 
dinal bending stresses thus obtaining 




We will then expect to obtain TOP> BOT and DITOP < DIBOT 
and we may already consider the respective expressions for 
DITOP, DIBOT and I with TS replaced by TOP 

Prom the diagram we may consider thefollowing equations: 
STH*DIBOT=SBH*DITOP 

and notice here that if we want to stress the material to the 
allowable design stresses we have 
STH = DST-STR3T 
SBH = DSC-STR3B 

we may combine these equations to obtain 

(DST-STR3T)*DIB0T = ( DSC-STR3B ) *DITOP 
STR3T and STR3B are given by 

S 2 

STR3T = 0.152222*HT0P*( ) c 

TOP 

S 2 

STR3B = 0.152222*HB0T*( ) 

BOT 
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anticipating here that S = Z*BOT because BOT is the thinner 
plating we have 




BOT 

and thus may write 

STR3B = 0. 152222*HB0T*Z 2 

Since 

S ? S 2 

( — r< ( — r 

TOP BOT 

because TOP > BOT, if we make 




TOP 

to simplify the equations that will follow we will understand 
that in doing so we are only increasing the value of STR3T. 
This as we see fron 

STH = DST-STR3T 

will make us have STH with a slightly lower value than the one 
that is allowed, so we are on the safe side when we consider 




TOP 



Now if we end up too much on the safe side and we have a TSTH 
too low we may bring it to a value closer to DST by considering 
a Z1 value as done and exemplified in the previous section and 
done again in this section for the case 1000 ft steel where 
for this particular exemple it is TSBH that needs correction 
as explained in detail . 



We may then write 



1 

Z = — ; - — 

•V I 1-P0I*P0I 

\ 0 . 30396 * *DSC 

1 E 

STR3T = 0. 152222*HT0P*Z 2 
STR3B = 0. 152222*HB0T*Z 2 
and from 

(DST-STR3T)*DIB0T = (DSC-STR3B)*DIT0P 
we have if we make 
A = DST-STR3T 
AB = DSC-STR3B 
P = (2*D+B)*T0P+B*B0T 

D*T0P* ( B+D ) D*(B*B0T+D*T0P) 

P P 

which gives after simplification 

A*<(B+D)-AB*D)*T0P = AB*B*B0T 
and making 

AC = ( A* ( B+D ) -AB*D ) / ( AB*B ) 
we have 

BOT = AC*T0P 

Following th e suggestions of a previous section we will 
also use the equation that will make the bottom to be stressed 
the most in compression under the larger of the moments being 
considered (HM in this case). 

We obtain 

SBH = HM*DIB0T/I = AB 

Replacing DIBOT and I by their expressions we obtain 
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HM*D*TOP*(B+D) 

AB o 5 

4*D^*T0P* ( ( 3*B*+2*B*D ) *BOT+ ( D 2 +2*B*D ) *TOP ) 

simplifying we end up with 

HM*(B+D) 

TOP 

AB*4*D*(D 2 +2*B*D) 

now making 

AD = HM*(B+D)/(AB*4*D*(D 2 +2*B*D) ) 

AE = (3*B 2 +2*B*D)/(D 2 +2*B*D) 
we have 

TOP = AD-AE*BOT 
substituting here 
BOT = AC*TOP 
we obtain finally 
AD 

TOP = 

1+AC*AE 

Considering here again the bottom as the location of the thinner 
plating (in this particular situation we are looking at), and 
knowing that it will be subject to a compressive stress of 
magnitude as large as DSC we may again write S = Z*BOT 
We finally obtain the total stress on top and bottom in both 
tension and compression by adding the longitudinal bending 
stresses to the tertiary stresses. 

TSTH = HM*DITOP/I+0.152222*HTOP*S 2 /TOP 2 
TSTS = SM*DITOP/I+0. 152222*HT0P*S 2 /T0P 2 
TSBH = HM*DIB0T/I+STR3B 

i 

TSBS = SM*DIB0T/I+STR3B 



3*B 2 +2*B*D 

— tj *BOT 

D +2*B*D 



NOTE: since usually the secondary stresses are assigned 
a constant value like 2000 psi or 3000 psi they can be a- 
ccounted for when giving a value to the design stresses. 
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As we see the results obtained for the 500 ft aluminum 
case are acceptable and agree with all our considerations in 
the derivation of the formulas, namely T0P>B0T. 

The following two runs, for 500 ft steel and 1000 ft alu- 
minum show that the results are in disagreement with T0P> B0T 
used in the formulas. We will see next why this happens and 
derive the proper set of formulas to be used for such cases. 
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T 0 P = 


? • 1023 


BOT = 


3 *8359 


s = 


13* .56*9 


DITOP= 


38 • * 76? 


D I BO T = 


3? • 95?* 


TSTHe 


25752.578 


TSTS = 


1^*27 .625 


STCRI T= 


19622.289 


T S p H = 


29999.980 


T S P S = 


2*5/48.805 


sbcri T = 


3 1’* 007 . 0 3 1 



• VlO OPERATING SYSTEM VERSION 1 REVISION 312 P3/t'*/7* GENERATED 



jOR HANDLING CHARGE 5 .35 / JOB .35 
30 LINES PRINTED PR? s 1.25 / v lN .15 
E8 CARDS READ * 1.5? / K CO *13 
P0 PLOTTER VECTORS S .2< s / j-D/ .07 
c 2 .'-ODEL 77 SECONDS s-25.70 / POUR .16 
00 NCDEl 87 SECONfS & 12.57 / Mr UR .07. 

TOTAL CHARGE * .79 



ElR A 9P 1*731 LOGGED OUT 75/01/7* 22:26. 



* 2 *. 80 LEFT AFTER 



7*/3fV7* 1 



03 LOGINS* 
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STEEL 

POI = 0.3 

E = 30*10 6 
AL = 500 
B a AL/8 =62.5 
DRA = B/3 = 20.833 
D = AL/9=55.556 

m = 1. 37l43*B*DRA*AL 2 = 44642. 1 87*1 0 4 
SM = B*DRA*AL 2 = 32551 . 562*1 0 4 
w = 0 . 30396 * (1-P0I 2 )*DSC/E = 276. 6*1 o " 6 
Z = 60.129 

HTOP = 8 
HBOT = 20.833 

STR3T = 0.152222*HTOP*Z 2 = 4402.865 note that this value 

should be corrected for the Z 
approximation and actually is 
slightly smaller 

STR3B = 0. 152222*HB0T*Z 2 = 11465. 6l 
We see from the tertiary stresses obtained that we end up with 
the following diagram of stress distributions: 
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From the diagram we see that theneutral axis will end up closer 
to the bottom and so we must consider BOT > TOP in the formulas. 



The same type of diagram will be obtained for 1000 1 aluminum case 
since the tertiary stresses obtained in such case are 
STR3T = 1498.73 
STR3B = 9873.078 

For the 1000' steel case the same type of diagram will occur 
due to the tertiary stresses being 
STR3T = 4402.86 
STR3B = 29004. 

Besides this situation requiring us also to consider the 
neutral axis closer to the bottom^ still another problem arises 
since STR3B is larger than DST=25000. This case is solved by 
properly correcting the inaccuracy that was introduced by the 
Z approximation as exemplified next. 
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Now let's derive the formulas capable of handling the 500' 
steel case, the 1000' steel case and the 1000' aluminum case. 

We must consider BOT > TOP and DIT0P>DIB0T 
DIT0P = D*(B+D)*B0T/P 
DIB0T = D*(B*T0P+D*B0T)/P 

I = 4*D 2 *B0T* ( ( 3*B 2 +2*B*D ) *T0P+ ( D 2 +2*B*D ) *B0T ) /P 

where 

P = 2*D*BOT+B*(TOP+BOT) 

Starting again with the equation 
STH*DIB0T = SBH^DITOP 
and making 

A = DST-STR3T 
AB = DSC-STR3B 
we obtain 

A*DIB0T = AB*DITOP 

Knowing that now the top is the thinner plating we anticipate 
S = Z*T0P 
and may write 

STR3T = 0. 152222*HT0P*Z 2 
STR3B = 0.152222*HBOT*Z 2 

S 

where the Z value in STR3B should be which has a smaller 

BOT 

value than Z. 

We then obtain: 

A*D*(B*T0P+D*B0T)/P = AB*D* ( B+D ) *B0T /P 



2l4 



Simplifying and making 

AC = (AB*(B+D)-A*D)/(A*B) 
we obtain 

TOP = AC'^BOT 

Making now the top to be stressed the most in tension we have. 
STH = HM*DITOP/l =A 

Replacing DITOP and I by their equations we have 

HM*D*(B+D)*BOT 

A 

4*D 2 *B0T*( ( 3*B b +2*B*1D)*T0P+(D 2 +2*B*D)*B0T ) 

and obtain 

HM*(B+D) 3*B 2 +2*B*D 

BOT *TOP 

A*4#D* ( D 2 +2*B*D ) D 2 +2*B*D 

Making 

AD = HM*(B+D)/(A*4*D*(D 2 +2*B*D) ) 

AE = (3*B 2 +2*B*D)/(D 2 +2*B*D) 
the formulas become 

BOT = AD-AE*TOP 
and substituting 

TOP = AC*BOT 
we finally have 
AD 

BOT = 

1+AC*AE 

The results obtained for the examples' chosen are presented next: 
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We could have obtained TSBH = 30000 by correcting the 
approximation used 2 = S/BOT as is done for the 1000' steel 
case that follows next. 

In this case we already anticipate the use of a Z1 
value as done in the previous section and according to the 
output obtained we will then obtain the correct result. 
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Now we see that with an assumed value Z1 = 35 we obtain 
S = 215.1808 

BOT = 3.686 so Z1 actual = 58.38 

TSBH = 4-7512. 793 (too high) 

TSBS also too high 

The Zl actual value obtained suggests that we should have 
assumed an higher Zl value and with an assumed value Zl = 4-0 
we then obtain 

s = 203.626 

BOT = 4-. 428 so Zl actual = 45.98 

TSBH = 3^131.855 (too high) 

TSBS also too high 

This suggests an assumed Zl value higher than 40 but now we may 
interpolate with the two results already obtained and determine 
the Zl value that will bring TSBH = DSC = 30000 psi and done 
properly this new Zl assumed will end up being Zl actual. 

From interpolation we obtain 

(40-35)^(34131. 855-30000) 

Zl = 40+ 41.544 

(47512.793-34131.855) 



and now assume this value obtaining the following result: 
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We see that TSBH is just a little too high and the same 
happens with TSBS; here it might be the case that TSBS should 
have been the one under consideration for the choice of a Z1 
value. Nevertheless for the purpose of this work this example 
shows the procedure to be used along with the equation derived. 
The accuracy obtained can be total. Here we assumed 
Z1 = 41.544 and ended up with Z1 actual = 41.83 which is close 
enough for our purpose and TSBH and TSBS is very close to what 
it should be to use the material effectively. 
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This last example shows that TSBH is low and the Z 
approximation used in the formulas here again should be 
corrected as already exemplified in order to bring it closer 



to the DSC value. 
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CONCLUSION 

Prom this work the author is led to conclude that there 
is no simplified approach to the neutral axis location in the 
midship section, capable of handling properly all or only 
a majority of cases. It is advisable to look into each parti- 
cular situation by first gaining some insight at what the 
diagrams of stress distributions may look like as exemplified 
in several sections. The two last sections exemplify well why 
this may be of some advantage and specially the last section which 
shows that while at first it was expected for the examples stu- 
died a neutral axis closer to the top, actually the results 
ended up requiring a neutral axis closer to the bottom and a 
new set of formulas had to be derived. 

As a future development for this work it is here sugges- 
ted that a more ship like cross section could be analysed as 
a step forward towards obtaining a method of determination of 
scantlings as direct as possible. This next step should consi- 
der the derivation of proper formulas for DITOP, DIBOT and I 
involving only the plate thicknesses as unknowns . It should 
be considered then, in order to simplify the problem. , a stan- 
dard type of frame and its scantlings should be constant throu- 
ghout the midship section. 
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